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Since  publication  of  Pourbaix's  classic  treatise  "Atlas  of  Electro¬ 
chemical  Equilibria  in  Aqueous  Solutions"  (Pergamon  Press,  Oxford, 
CEBELCOR ,  Brussels,  1966),  there  has  been  interest  in  computerizing  the 
calculation  of  Pourbaix  diagrams  in  order  to  avoid  the  tedious  task  of 
calculating  these  diagrams  "long-hand"  for  complex  systems.  Reference  1 
reports  a  computerized  method  developed  at  University  of  Florida  and  cur¬ 
rently  in  use  in  a  number  of  laboratories  which  have  access  to  main-frame 
computing  facilities.  On  the  assumption  that  personal  computers  can  be 
made  more  readily  available  to  scientists  and  engineers,  there  was  inter¬ 
est  in  adopting  the  existing  program  for  use  on  an  Apple  XI  personal 
computer . 

ONR  Contract  #N00014-82-C0718  provided  support  for  this  effort.  The 
new  program  developed  at  University  of  Florida  has  the  following  capabil¬ 
ities  and  advantages: 

(1)  The  program  is  written  in  UCSD-Apple  Pascal  which  makes  it  adaptable 

to  many  microcomputer  systems. 

(2)  The  method  used  to  balance  the  equilibrium  chemical  equations  allows 


the  consideration  of  complicated  metal-ion-water  systems. 

(3)  The  program  can  be  easily  modified  for  high  temperature  applications' 

through  the  inclusion  of  an  available  subprogram.  £ 

(4)  The  method  used  to  determine  areas  of  predominance  is  intuitively  ^ 

obvious,  rigorously  correct,  and  continuum  in  its  approach.  n 

(5)  The  program  can  be  applied  to  determine  equilibrium  potential-pH  Ajf 


diagrams  for  even  very  complex  syst 
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The  Fe-Cl-HjO  system  provides  a  simple  example  of  the  application  of 
the  program.  The  metallic  species  which  will  be  considered  are  shown  in 
Table  1 .  Later,  a  more  complex  example  will  be  given. 


Table  1  .  DATA  EMPLOYED  IN  COMPUTERIZED  CALCULATION  OP  THE  POTENTIAL  vs 
pH  DIAGRAM  FOR  THE  Pe-Cl-HjO  SYSTEM  AT  25°C 


Species  considered 
Chemical  Computer 

notation  notation 

Gibb's  free  energy 
of  formation 
(in  cals/M) 

Assumed 

thermodynamic 

activity 

Cl- 

CL- 

-31372.0 

1.0  x  10-1 

Fe 

FE 

0.0 

1 .0 

Fe3D4 

FE304 

-242700.0 

1.0 

Fe2°3 

FE203 

-177400.0 

1 .0 

HFeOj" 

HFE02- 

-90300.0 

1.0  x  10“6 

FeOH2+ 

FEOH++ 

-54830.0 

1.0  x  10"6 

Fe(OH)2" 

FE(0H)2+ 

-104700.0 

1.0  x  10~6 

FeCl2(Aq) 

FBCL2(AQ) 

-81590.0 

1 .0  x  10“6 

FeCl2+ 

FBCL++ 

-34400.0 

1.0  x  10-6 

Since  computer  line  printers  display  only  capital  letters,  it  is 
necessary  to  develop  a  suitable  convention  for  writing  chemical  equations 
entirely  on  one  line  (no  subscripts )j and  in  capital  letters.  Any  reason¬ 
able,  internally  consistent,  convention  nay  be  employed.  Further 
reference  to  Table  1  will  supply  insight  into  the  type  of  convention  em¬ 
ployed  by  the  authors  in  this  example.  In  order  to  avoid  confusion  in 
reading  the  balanced  chemical  equations,  the  computer  places  an  asterisk 
between  all  coefficients  and  the  chemical  species  to  which  they  refer. 
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Input  of  Data 


Data  is  entered  into  the  microcomputer  interactively.  The  user  is 
prompted  for  each  piece  of  data  to  be  entered.  In  order  to  prevent 
crashing,  all  data  is  read  as  characters  and  converted  to  integer  or  real 
values  as  required. 

The  first  piece  of  data  to  be  entered  is  the  temperature,  in  degrees 
Kelvin,  to  be  used  in  Nernst  equation  calculations.  Temperature  correc¬ 
tions  of  free  energy  data  may  be  done  using  a  separate  subprogram  based 
on  the  methods  of  Criss  £  Cobble,  (ref.  2) 

The  number  of  non-metallic  ionic  species  to  be  considered  is 
entered.  In  this  example  the  number  1  is  entered,  since  only  Cl~  is  con¬ 
sidered.  Pertinent  ionic-species  data  in  the  form  of  each  species'  name, 
free  energy  of  formation,  activity,  electrical  charge  and  elemental  com¬ 
position  is  then  input. 

.  Immediately  after  the  ionic  species  data  the  number  of  metallic 
species  to  be  considered  is  entered  (for  the  Fe-Cl-HjO  system  this  number 
is  8,  since  Cl~  mould  not  be  counted  as  a  metallic  species).  Other  per¬ 
tinent  data  similar  to  that  listed  above  for  the  ionic  species  is  then 
input  for  each  metallic  species. 

When  data  entry  is  completed  the  user  is  asked  if  the  data  requires 
correction.  The  user  is  able  to  select  and  modify  any  piece  of  data. 

The  data  is  written  to  the  terminal  or  printer  as  specified  by  the 
user.  The  tabular  fora  allows  checking  for  accuracy  and  the  user  is 
again  offered  the  opportunity  to  modify  the  entered  data. 
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Chemical  equation  balancing  and  Equilibrium  Line  Determination 

When  considering  the  possible  formation  of  N  different  metallic 
species 

N  (N  -  1) 

2 

chemical  equilibrium  equations  exist.  The  program  balances  chemical 
equations  between  all  of  the  possible  pairs  of  metallic  species.  The 
program  takes  all  possible  pairs  of  metallic  species  and  balances  the 
equations  in  terms  of  water,  hydrogen  ions,  electrons,  etc.  The  balanced 
chemical  equations  are  used  by  the  computer  to  determine  and  tabulate  the 
equilibrium  line  equations  by  means  of  Nernst  equation  calculations. 

(For  a  simple  explanation  of  the  calculations  Involved,  consult  Ref.  3) 

Areas  of  Predominance 

Determining  the  areas  of  predominance  on  a  potential-pH  diagram  re¬ 
quires  a  proper  understanding  of  the  information  conveyed  by  each  line  on 
the  diagram.  Each  line  represents  equilibrium  coexistence  between  two 
metallic  species  (with  certain  assumptions  regarding  ionic  activities). 
When  any  line  is  crossed  one  metallic  species  becomes  more  energetically 
favorable  than  another.  Each  line,  therefore,  delineates  a  boundary  be¬ 
tween  two  adjacent  areas  on  the  diagram.  Each  area  represents  the  region 
(in  potential  vs  pH  coordinates)  in  which  a  particular  species  is  thermo¬ 
dynamically  stable  and  dominant.  By  simultaneously  considering  all  of 
the  lines  which  correspond  to  equilibria  involving  a  given  species,  an 
area  is  traced  out  in  which  the  species  in  question  is  "not  excluded". 
This  then  is  the  area  of  predominance  for  that  species. 
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When  each  species  is  considered  (one  at  a  tine)  in  the  above  nanner, 
a  series  of  separate  areas  of  predoninance  results.  When  these  separate 
areas  are  fitted  together,  they  forn  the  finished  theoretical  potential- 
pH  diagram  without  any  overlaps. 

A  Mathematical  analog  can  be  drawn  by  considering  the  equilibrium 
lines  as  inequalities.  For  example,  in  considering  the  equilibrium 
between  Fe  and  Fe3<>4,  Fe  is  the  predominant  species  at  values  of  poten¬ 
tial  and  pH  less  than  those  described  by  the  equilibrium  coexistence  line 
between  Fe  and  Fe^O^  (hence  the  term  "inequality").  Thus  a  given  species 
will  exist  whenever  the  potential  and  pH  satisfy  all  of  the  line  inequal¬ 
ities  pertaining  to  the  species.  If  no  potential  and  pH  simultaneously 
satisfies  all  of  the  inequalities  of  a  given  species,  the  species  has  no 
area  of  predominance. 

The  field  of  operations  research  has  been  developed  by  industrial 
and  systems  engineers  to  deal  with  such  mathematical  analogs..  In  opera¬ 
tions  research  the  area  which  mutually  satisfies  a  set  of  linear  inequal¬ 
ities  is  termed  a  "convex  polygon".*  Methods  of  determining  convex 
polygons,  therefore,  lend  themselves  to  determining  areas  of  predoninance 
on  potential-pH  diagrams. 

Computer  determination  of  areas  of  predominance 

The  equation  balancing  portion  of  the  program  is  designed  to  write 
all  electron-transfer  equations  as  oxidation  reactions.  This  results  in 


*  Convex  polygon -an  area  such  that  if  a  line  segment  is  drawn  between  any 
two  arbitrarily  selected  points  in  the  polygon  all  points  on  the  line 
segment  will  also  lie  within  the  polygon.  (Kef.  3) 
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two  important  consequences.  First,  for  all  constant-pH  lines  the 
metallic  species  on  the  lower  pH  {left)  side  of  the  line  will  appear  on 
the  left  side  of  the  chemical  equation.  Second,  for  all  non-vertical 
lines  the  metallic  species  which  lie  on  the  lower  potential  sides  of  the 
lines  will  also  be  on  the  left  side  of  the  chemical  equation.  By  noting 
whether  the  various  lines  are  vertical  or  nonvertical  and  on  which  side 
of  the  chemical  equations  the  species  occur,  lines  can  be  sorted  into  the 
following  four  categories:  (1)  minimum  pH  lines;  (2)  maximum  pH  lines; 
(3)  potential,  upper  bound  lines;  and  (4)  potential,  lower  bound  lines. 
The  program  uses  this  sorting  technique  to  determine  the  convex  polygon 
for  each  metallic  species  which  simultaneously  satisfies  all  the 
equilibrium  line  equations  of  each  individual  species.  The  resulting 
polygons  are  displayed  to  the  user  by  the  computer  as  a  set  of  points 
which  constitute  the  corners  of  each  convex  polygon.  Plotting  the  deter¬ 
mined  points  (expressed  in  potential  and  pH  coordinates)  and  connecting 
them  with  straight  lines  results  in  a  diagram  like  the  one  shown  in  Fig. 

1  for  Fe-Cl-HjO. 


FIG.  1.  Potential  vs  pH  diagram  for  the  Fe-CI-H_0  system  at  25*C 

as  determined  by  computer  considering  tne  species  shown  in 

Table  1 .  /■ 
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Determination  of  more  Complex  Systems 


More  complex  metal-ion-water  systems  involving  ionic  species  which 
Calculation  of  potential-pH  diagrams  of  metal-ion-water  systems  by  computer 


PIG.  2.  Potential  vs  pH  diagram  for  the  Fe-Cl-S-HjO  system  at 

2S*C  as  determined  by  the  computer  considering  the  species 
shown  in  Tables  1  and  2. 

•  2— 

have  their  own  potential-pH  equilibria  (e.g.  HjS,  HS  ,  S  )  appear  to 
violate  the  previous  theory  for  determining  regions  of  predominance. 

Such  is  not  the  case  when  the  program  is  properly  employed. 

Addition  of  the  species  shown  in  Table  2  to  those  in  Table  1 
results  in  the  computer  determined  diagram  shown  in  Pig.  2.  Note  even 
though  the  region  of  predominance  for  PeS  does  not,  as  a  whole,  consti¬ 
tute  a  convex  polygon,  the  portions  of  the  region  which  lie  in  each  of 
the  individual  areas  of  predominance  for  HjS,  HS-  and  S2~  respectively, 
are  in  each  case  convex.  Thus  when  it  becomes  necessary  to  formulate 
such  e  diagram  it  requires  the  initial  computer  calculation  of  sore  than 


one  diagram.  In  this  case  three  diagrams  were  determined  considering 


—  2— 

individually  H^S,  HS  and  S  The  appropriate  portion  of  each  of  the 
three  diagrams  was  then  used  to  construct  the  final  diagram 


Table  2.  DATA  EMPLOYED  IN  COMPUTERIZED  CALCULATION  OF  THE  POTENTIAL  vs 
pH  DIAGRAM  FOR  THE  Fe-Cl-S  SYSTEM  AT  25°C 


Species  considered 

Chemical  Computer 

notation  notation 

Gibb's  free  energy 
of  formation 
(in  cals/M) 

Assumed 

thermodynamic 

activity 

h2s 

H2S 

-6660.0 

e 

o 

X 

— * 

o 

1 

HS" 

HS- 

2880.0 

1  .0  x  10-1 

s2- 

ss — 

20500.0 

T 

o 

X 

o 

t 

FeS  (Solid) 

FES  (SOL) 

-24000.0 

1  .0 

FeS2  (Solid) 

FES2  (SOL) 

-39900.0 

1  .0 

FeSQ4  (Aq) 

FES 04  (AQ) 

-196820.0 

1.0  x  10"6 

shown  in  Fig.  2  taking  into  account  the  domain  of  stability  of  each 
ionic  species.  (For  simplicity  the  diagrams  are  drawn  with  straight- 
line  segments  meeting  at  the  "equal  dominance"  lines  for  adjacent  ionic 
species.  Actually  at  such  lines  (e.g.  KS~/S2~)  the  activities  of  both 
species  are  equal  at  the  equal  dominance  line  but  one  species  rapidly 
predominates  as  one  moves  away  from  the  line.)  Using  this  technique 
even  more  complicated  systems  can  be  determined. 
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Summary  of  Program  Capabilities 


The  interactive  microcomputer  version  of  the  pH  potential  program 
is  largely  a  translation  of  the  WATFIV  version  and  takes  advantage  of 
UCSD  Pascal  features. 

The  described  computer  program  is  arbitarily  designed  to  consider  a 
maximum  of  35  metallic  species  and  10  non-metallic  ionic  species  simul¬ 
taneously.  Redimensioning  of  the  program  matrices  will  allow  for  any 
desired  increase  or  decrease  in  the  program's  capacity. 

The  structure  of  the  program  and  the  UCSD  operating  system  facili¬ 
tate  the  addition  of  extensions  to  the  program  through  the  independent 
development  of  new  modules.  It  is  hoped  that  support  can  be  obtained 
for  development  of  additional  modules  including  a  data  maintenance  rou¬ 
tine  that  will  record  on  disk  a  library  of  data  for  selected  metal-ion- 
water  systems  and  an  accompanying  graphics  package  that  will  permit 
interactive  terminal  and  printer  output  of  generated  pH  potential 
diagrams . 
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Attachments  Include : 


Flow  chart 

Modular  Documentation 

The  Main  Program  with  Accompanying  Units 
User's  Guide 
Example  of  Output 
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PH  -  POTENTIAL  PROGRAM 


MODULAR  DOCUMENTATION 


MODULAR  DESCRIPTION  FORM 


A.  MODULAR  NAME  -  Metal  input  (*  part  of  unit  in  data*) 


B.  MODULE  FUNCTION  -  Prompts  user  for  metallic  species  data 


C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Terminal: Number  of 
metallic  species  (NP) 

Name  of  metallic  specie  (SP  (M,  HH] ) 

Free  energy  of  formation  (F  [M] ) 

Ionic  activity  (AM  [M] ) 

Electrical  charge  (Z  [Mj ) 

Number  of  metal  atoms  (Metal  (M) ) 

Number  of  oxygen  atoms  (Oxygen  (M] ) 

Number  of  hydrogen  atoms  (Hydro  (Ml) 

Number  of  characteristic  atoms  for  each  ionic  species  (ION  [I,  Ml ) 


D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  None 


E.  SUMMARY  OF  PROCESS  -  Sequentially  prompts  user  for  metallic  species 
data.  All  real  values  read  as  characters  by  the  procedure  inreal, 
and  are  converted  to  real  by  CHAR  to  REAL.  Integer  values  are  also 
read  as  characters  and  are  converted  to  integer  by  CHAR  to  I NT. 

Name  of  species  is  read  in  by  the  procedure  GETMET. 


F.  CALLS  THESE  MODULES  -  Procedures: INCHAR,  INREAL,  GETMET,  LF. 
Unit: Convert 


G.  CALLED  BY  THESE  MODULES  -  Main 


MODULAR  DESCRIPTION  FORM 


A.  MODULAR  NAME  -  Data  input  (*  part  of  unit  in  data  *) 


B.  MODULE  FUNCTION  -  Prompts  user  for  non-metal lie  ionic  species  data 


C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Terminal : Temperature 
(T) 

Free  energy  of  formation  of  HjO  (FEHjO) 

Number  of  non-metallic  ions  (NION) 

Name  of  species  (Name  [K,  HHj ) 

Free  energy  of  formation  (FION  (1)) 

Ionic  activity  (AION  [1]) 

Electrical  charge  ( ZION  [ 1 ] ) 

Number  of  characteristic  atoms  (BION  [1)) 

Number  of  Oxygen  atoms  (OXY  [1]) 

Number  of  hydrogen  atoms  (HYD  (1]) 


D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  None 


E.  SUMMARY  OF  PROCESS  -  Sequentially  prompts  user  for  non-metallic 

ionic  species  data.  All  real  values  are  read  as  characters  by  the 
procedure  inreal  and  converted  to  real  by  CHAR  to  REAL.  Integer 
values  are  also  read  as  characters  and  are  converted  to  Integers  by 
.CHAR  to  I NT.  Name  of  species  is  read  by  the  procedure  GETNONMET. 


F.  CALLS  THESE  MODULES  -  Procedures t INCHAR ,  INREAL,  GETNONMET,  LF. 
Units Convert. 


G.  CALLED  BY  THESE  MODULES  -  Main 


MODULAR  DESCRIPTION  FORM 

A.  MODULAR  NAME  -  Initialize 

B.  MODULE  FUNCTION  -  Initialization  of  index  values,  ionic  species 
values,  and  species  name  array 


C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  None 

D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  None 

E.  SUMMARY  OF  PROCESS  -  Initializes  program  counter  index  to  1 . 
Initializes  FION  [1],  OXY  [1],  to  0.  Initializes  S,  Name,  X,  SP  to 
all  blanks  ("b"). 

F.  CALLS  THESE  MODULES  -  None 

G.  CALLED  BY  THESE  MODULES  -  Main 


MODULAR  DESCRIPTION  FORM 


A.  MODULAR  NAME  -  Main  program 


B.  MODULE  FUNCTION  -  Controls  execution  of  program 


C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Interactively  prompts 
user  for  display  option  (printer  or  terminal) 


D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  None 


E.  SUMMARY  OF  PROCESS  -  Sequentially  calls  primary  program  modules  for: 

1 .  input  of  data 

2.  correction  of  data 

3.  display  of  data 

4.  correction  of  data 

5.  equation  balancing 

F.  CALLS  THESE  MODULES  -  Initialise,  data  input,  metal  input,  nonmetal 
display,  metal  display,  data  correct,  Bg  balancer 


G.  CALLED  BY  THESE  MODULES  -  None 
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MODULAR  DESCRIPTION  FORM 


A.  MOOULAR  NAME  -  Stable  1  (*  part  of  unit  stable  * ) 


B.  MODULE  FUNCTION  —  Determining  areas  of  predominance 


C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  &  TERMINAL)  -  None 

D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  &  TERMINAL)  -  Listing  of  areas  of 
predominance,  line  segment  data.  Listing  of  minimum  and  maximum 
constant  pH  line  data. 


E.  SUMMARY  OF  PROCESS  -  Determine  minimum  pH  line  for  metallic  species 
(Stable  3).  Determines  maximum  constant  pH  line  for  metallic 
species  (Stable  4).  Calculates  upper  bound  lines  to  determine  the 
line  which  supplies  minimum  potential  for  upper  bound  of  stability 
region  at  minimum  pH  (UPBNDLNS).  Calculates  intersection  points 
between  previously  determined  upper  bound  lines  and  remaining  upper 
bound  lines  (INTPTSHI) .  Calculates  lower  bound  lines  (LOWBNDLNS) 
and  intersection  points  for  lower  bound  lines  (INTPTSLO).  Tests 
upper  and  lower  bounds  to  determine  whether  they  intersect  within 
the  minimum  and  maximum  pH  lines.  If  they  intersect,  the  points  are 
determined  and  data  is  printed.  Orders  pH  values  at  which  either 
two  upper  bound  lines  or  two  lower  bound  lines  intersect  from  PHMIN 
to  PHMAX.  Finds  point  where  upper  bound  exceeds  lower  bound. 
Determines  intersection  point  on  maximum  pH  side.  Prints  data. 


F.  CALLS  THESE  MODULES  -  Procedures : Start ,  Set  flags.  Stable  2,  Stable 
3,  stable  4,  UPBNDLNS,  INTPTS  HI,  TEMP  1,  LOWBNDLNS,  INTPTSLO,  TEMP 
2,  Stable  5,  Stable  6,  Stable  7,  stable  8,  Stable  9,  Stable  10, 
Stable  1 1 . 


G.  CALLED  BY  THESE  MODULES  -  Eq  Balancer 


MODULAR  DESCRIPTION  FORM 


A.  MODULAR  NAME  -  Eq  Balancer 


B.  MODULE  FUNCTION  -  Balances  equations  between  all  possible  pairs  of 
metallic  species. 


C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  None 


D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Terminal  or 
Printer:  Balanced  equations.  Equilibrium  line  equation. 


E.  SUMMARY  OF  PROCESS  -  Balances  chemical  equations  between  all 

possible  pairs  of  metallic  species.  Takes  all  pairs  of  metallic 
species  and  balances  them  in  terms  of  water,  hydrogen  ions, 
electrons,  etc.  NERNST  equation  calculations  are  then  used  to 
tabulate  the  equilibrium  line  equations.  Procedure  diagram  then 
calls  stable  1 . 


F.  CALLS  THESE  MODULES  -  Procedures : EQB1  ,  SUMNUM,  Update,  LOGSUM, 
FINDSIGN,  EQB4,  CALCEQEQ, 

Diagram,  Stable  1 . 

G.  CALLED  BY  THESE  MODULES  -  Main 


MODULAR  DESCRIPTION  FORM 


A.  MODULAR  NAME  -  Nonmetcor  (*  part  of  unit  in  data  *) 


B.  MODULE  FUNCTION  -  Correction  of  non -metal lie  species  data  entered. 


C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  fi  TERMINAL)  -  Terminal:  Parameter 
to  be  modified.  (Name,  FION,  AION,  ZION,  BION,  OXY,  HYD) .  New 
value  for  parameter  specified. 


D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  S  TERMINAL)  -  Terminal:  Menu  of 
parameters.  Old  value  of  parameter. 


E.  SUMMARY  OF  PROCESS  -  Asks  which  ion  to  modify.  Displays  menu  of 
options  to  be  modified. 

1 .  ion  name 

2.  free  energy  of  formation 

3.  activity 

4 .  charge 

5.  number  of  characteristic  atoms 

6.  number  of  oxygen  atoms 

7.  number  of  hydrogen  atoms 

Writes  old  value  of  parameter,  prompts  for  and  reads  new  value. 


F .CALLS  THESE  MODULES  -  Procedures :NMETC1 ,  NMETC2,  NMETC3,  NMETC4, 
NMETCS,  M1ETC6,  NMETC7 , 

6ETN0NMET,  INREAL,  INCHAR. 

UnitiConvert  . 


G. CALLED  BY  THESE  MODULES  -  Data  correct 


MODULAR  DESCRIPTION  FORM 


A.  MODULAR  NAME  -  Meteor 


B.  MODULE  FUNCTION  -  Correction  of  metallic  species  data  entered. 


C.  INPUTS  ( ARGUMENTS ,  COMMON,  FILE  £  TERMINAL)  -  Terminal :speci as  and 
parameter  to  be  modified  (SP,  F,  AM,  z,  Metal,  Oxygen,  Hydro,  Ion). 

New  value  for  parameter  specified. 


D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Terminal : Menu  of 
parameters.  Old  value  of  parameter . 


E.  SUMMARY  OF  PROCESS  -  Asks  for  species  to  be  modified.  Displays  menu 
of  options. 

1 .  specie  name 

2.  free  energy  of  formation 

3.  activity 

4.  charge 

5.  number  of  metal  atoms 

6.  number  of  oxygen  atoms 

7.  number  of  hydrogen  atoms 

8.  number  of  characteristic  atoms 

Writes  old  value  of  parameter,  prompts  for  and  reads  new  value. 


F. CALLS  THESE  MODULES  -  Procedures « METC1 ,  METC2,  METC3 ,  METC4,  METC5 , 
METC6,  METC7,  METC8, 

GETMET,  INCHAR,  INREAL. 

Unitt Convert 


G. CALLED  BY  THESE  MODULES  -  Data  correct 


MODULAR  DESCRIPTION  FORM 

A.  MODULAR  NAME  -  Fecor  (*  part  of  unit  in  data  *) 


B.  MODULE  FUNCTION  -  Correction  of  free  energy  of  formation  of  H20 
entered. 


C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Terminal:  FEH20 


D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Old  FEH20  value 


E.  SUMMARY  OF  PROCESS  -  Prints  current  value  of  FEH20.  Prompts  for  and 
reads  new  value  of  FEH20. 


F.  CALLS  THESE  MODULES  -  Procedure: Inreal 
Unit:Convert 


G.  CALLED  BY  THESE  MODULES  -  Data  correct 


MODULAR  DESCRIPTION  FORM 


A.  MODULAR  NAME  -  Tempcor 

B.  MODULE  FUNCTION  -  Correction  of  temperature  entered 

C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Terminal:  T 

D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Old  T  value 

E.  SUMMARY  OF  PROCESS  -  Prints  out  current  value  of  temperature. 
Prompts  for  and  reads  new  value  of  temperature. 

F.  CALLS  THESE  MODULES  -  Procedure: Inreal 
Unit: Convert 

G.  CALLED  BY  THESE  MODULES  -  Data  correct 


MODULAR  DESCRIPTION  FORM 


A.  MODULAR  NAME  -  Data  correct  (*  part  of  unit  in  data  *) 


B.  MODULE  FUNCTION  -  User  correction  of  data  entered 


C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Data  type  to  be 
modified  (Opt.) 


D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  None 


E.  SUMMARY  OF  PROCESS  -  Prompts  user  for  type  of  data  that  requires 
modification.  The  options  are: 

1 .  temperature 

2.  free  energy  of  t^O 

3.  non -metallic  species  data 

4.  metallic  species  data 

5.  data  is  now  correct 
Calls  appropriate  procedure. 


F.  CALLS  THESE  MODULES  -  Procedures:  Tempcor,  Fecor,  Nonmetcor,  Meteor 

G.  CALLED  BY  THESE  MODULES  -  Main 


MODULAR  DESCRIPTION  FORM 


A.  MODULAR  NAME  -  Metal  display 

B.  MODULE  FUNCTION  -  Displays  metallic  species  data 

C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  6  TERMINAL)  -  None 

D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  £  TERMINAL)  -  Terminal  or 

Printer:  SP,  F,  AM,  Z,  Metal,  Oxygen,  Hydro,  Ion 

E.  SUMMARY  OF  PROCESS  -  Displays  metallic  species  data  in  tabular  form 
to  terminal  or  printer  as  specified  by  user. 

F.  CALLS  THESE  MODULES  -  Procedures:  Dis  1  MET,  Dis  2  MET,  Dis  3  MET, 
Draw. 

G.  CALLED  BY  THESE  MODULES  -  Main 


MODULAR  DESCRIPTION  FORM 

A.  MODULAR  NAME  -  Nonmetal  display 

B.  MODULE  FUNCTION  -  Displays  non-metallic  species  data 

C.  INPUTS  (ARGUMENTS,  COMMON,  FILE  6  TERMINAL)  -  None 

D.  OUTPUTS  (ARGUMENTS,  COMMON,  FILE  &  TERMINAL)  -  Terminal  or 

Printer:  T,  FEH20,  Name,  FION,  AION,  ZION,  BION  OXY,  HYD 

E.  SUMMARY  OF  PROCESS  -  Displays  non-metallic  species  data  in  tabular 
form  to  terminal  or  printer  as  specified  by  user. 

F.  CALLS  THESE  MODULES  -  Procedures:  ions  out  1,  ions  out  2,  draw. 

G.  CALLED  BY  THESE  MODULES  -  Main 


1 


<  *  SS**  ) 

PROCRAM  PH POTENT I AL  <  I NPUT , OUTPUT >  ; 

TYPE  TRANSCEND ' CETVALUE , CONVERT , INDATA , STABLE ; 

AONAME  »ARHAV[ 1.15,1  15]  OF  CHAR, 

MATR I  X 8 0  =  ARRAY C  1  .  .801  OF  CHAR; 

AO  I NT 1 0  s ARRAY  Cl.  .103  OF  INTEGER; 

?£5.2,X5V5*ARR*YC  ‘  23  0r  ARRAYC1  595  3  OF  REAL 

AOCH15  = ARRAY Ct  153  OF  CHAR; 


VAR 

C , COE  F 

JINDEX , INDEX 
MUN , A , NPP 
NUMB ,NUM ,NUME 
I  I  , HH , KK 
J,K,L, I ,M 
IONS , NON I  ON 
SAVE , TEMP , S 
MAT 
X 

NUMF , SUMAIO 
ACT, ACT1 , ACT2 
DELC , PH 
HE , ZAPl , ZAP2 
EN, ZAP3 , ZAP4 
OPT1 , NOPT 
AB 


AO 115; 
INTEGER; 
INTEGER ; 
INTEGER ; 
INTEGER ; 
INTEGER ; 
INTEGER ; 
AOCH1 5 ; 
MATR 1X80 ; 
AONAME ; 
REAL  ; 

REAL  ; 

REAL  ; 

REAL  ; 

REAL  , 

CHAR  ; 
BOOLEAN; 


PROCEDURE  INITIALIZE; 

BEGIN 

INDEX  : =1  ; 

FIONC 1 3  : «0  0  ; 

AIONtll ;«1  0; 

ZIONC 1 3  : =0  ; 

BIONC 1 3  : =0  ; 

OX  Y  C 1  3  : =0  ; 

HYDC 13*0; 

FOR  I : =1  TO  15  DO 
BEGIN 

SC  I  3  :  =  •  •  ; 

FOR  J : *  1  TO  15  DO 
BEGIN 

NAME C I , J3 : - •  • ; 

Xt I  , J3  : -  1  •  ; 

END; 

END; 

FOR  I : -1  TO  35  DO 
FOR  J:«l  TO  15  DO 
SPC I , J] : - •  • ; 

END; 


PROCEDURE  FINDS1GN; 

BECIN 

FOR  J : -1  TO  NUMB  DO 
IF  C0EFCJ3OD  THEN 
BEGIN 
I  :  -  I  1  ; 

S  t  I  3  .  .  •  ♦  ’  ; 

IF  1-1  THEN  SC  I  3  : •  ; 

IF  COEFCJ3<0  THEN  S  C I  3  : _ •  - •  - 
IF  J-CNION+3)  THEN 
BECIN 

SC  I  3  : - • « •  ; 

CC  I  3  -COEFC J3  ; 

END 

ELSE  CC 13 : -ABSC COEFC J3 > ; 

FOR  HH : ■ 1  TO  IS  DO  X C I , HH 3 ; -NAME C J , HH 3 ; 


END  ; 


PROCEDURE  UPDATE; 

BEGIN 

FOR  I : = 1  TO  NUMB  DO 
COE  F  C  I  3  : --COEFC I  3  . 

J : -COEFC 1 3  ; 

FOR  HH:=1  TO  15  DO 
BEGIN 

NAME  £ 1 ,HH3 : -NAME  INI  ON* 3 ,  HH3  , 
NAME  C  N I ON  +  3 , HH3  :  -TEMP  C  HH 3  ; 
END  ; 

COE Ft  1 3  : --COEFCNI ON  +33  ; 

COE  F  C  NI ON  + 3  3  : --J; 

ACT ; -ACT 1 ; 

ACT 1 : -ACT2 ; 

ACT2 : -ACT; 

LINEC 1 ,MI 3 ; -A; 

LINEC2 , M I  3  :  =M ; 

DE  LG : --DELC ; 

END  ; 

PROCEDURE  SUMNUM; 

BEGIN 


MUN: -0 ; 

NUM: =0 ; 

NUMF : -0 ; 

FOR  I : -1 

TO  NION  DO 

BECIN 

NUM: 

-NUM  +  COEF C I  +  2  3  *OXYC I 3 

MUN : -MUN  +  COEF  C I +  2 3  *KYDC I  3  ; 

NUMF  ;  -NUMF + COEF  CI+23*FI  ONC  I  3  ; 

END; 

COEFC  23  ; -OX YGEN C A 3 «METAL CM 3 -OX YGENCM3 «METAL C A 3 -NUM ; 

COEFCNION+43  :  -COEF  £23*2  -HYDRO C  A  3  *COEF  CNI ON+  3  3  +HYDROCM3  *COEF  t  1  3  +MUN 
DELG:  -COEF  C  NI ON*  3  3*FCA3-COEFt 1 3  * F C M3 -COEF C 2 3 « FEH20-NUMF ; 

NUME : -0 ; 

FOR  I  :  -1  TO  NION  DO 

NUME : -NUME ♦ COE F 1 1*23*21 ONC I  3  ; 

COEFCNION+53  ; -COEFCNION  +  4  3 -NUME -COEFC 1 3  « Z CM) +COEF CNI ON+ 3 3  « Z C A3  ; 
END; 

PROCEDURE  EQWRITE; 

BEGIN 
KK  •  -  1  ; 

WHILE  MATCKK3 O'* 1  DO 
BEGIN 

VR I TE ( F  F . MAT  C  KK  3 )  ; 

KK : -KK+1 ; 

IF  MATCKK3- *  -  *  THEN 
BEGIN 

VRITELNC  FF .MATCKK3 )  ; 

KK : -KK+ 1 ; 

WRITE ( FF , •  • 

WRITE<  FF ,  '  ' ) ; 

END; 

END; 

LNFD  < 1 )  ; 

WRITELNI FF , *  ! 

END; 


PROCEDURE  SUB 1 0 ; 

BEG  IN 

LTYPECMI 3 : -0 ; 

EQUATC 1  .MI  3  : -0 ; 

EOUATC  2 .MI  3 : -0 ; 

WRITE ( FF , '  '.MI: 3,'  !  NO  LINE  GENERATED'); 

WRITE  <  FF ,  '  !  ’>; 

EQWRITE ; 

END; 

PROCEDURE  SUB  1 1  ; 

BEGIN 

ZAP  1  ; -DELG/ (2 . 30  3*1  . f«7«T) ; 

ZAP2 : -COEF  CNI ON+3  3  *  LOG ( ACT 2 > ; 

ZAP3 ; -COEFC 1 3  *LOC( ACT1 >  ; 

PH:«(ZAP1 ♦ Z AP 2  - Z  A P 3 - BUM A  10) /COEFCNION+43 ; 
LTYPECMI 3  -l ; 

EOUATC 1 ,MI 3 : -PH; 

EQUATC  2 . MI  3 : -0  ; 

WRITE  <  FF ,  '  '  , MI  : 3 ,  '  !  PH-', PH  1:2); 

WR ITE  <  FF ,  '  !  '); 

EQWRITE ; 

END; 


I  '  )  ; 


!  '  )  ; 


PROCEDURE  CALCEQEQ, 

BEGIN 

I F < COEF CNION+5 I <0 )  AND  (COEFCNI ON 443-0)  THEN  SUB10 
ELSE 

IF  COEFCNION+5 I =0  THEN  SUBI1 
ELSE 
BECIN 

Z  A  P 1  :  = I  . ?8 4 4  *T/ COEFCNI ON* 5  1 ; 

Z AP2 : >DELG / (COEFCNI ON* 3 1  *  Z  30  6  0  .  0)  ; 

ZAP3  : -COEFCNI ON*  3 1  *  LOG  <  ACT2 >  ; 

ZAP 4  : «  COE  F  C 1 3  *LOG  <  ACT1 )  ; 

HE:»-< COEFCNI ON+41«ZAP1*0. 0001); 

EN: =ZAP2  +  ZAP1 * (ZAP3-ZAP4 -SUM A I O ) *  0  .0001  ; 
LTYPECMI 3=2; 

EQUATC 1  ,MI  3  ; =EN ; 

EOUATC 2 , MI  3  : -HE  ; 

IF  HE  =  0  THEN 
BEGIN 

WR I TE  <  FF ,  *  •.MI^,1  •  E-’,EN;8;4); 

WRITECFF , ‘  !  • > ; 

EQWRITE ; 

END 

ELif  HE<0  THEN 
BECIN 

HE  ;  = -HE ; 

WRITE  <  FF  ,  '  !  E=  '  ,  EN  :  8  :  4  )  ; 

WRITE  <  FF ,  *  -  *  , HE : 7  : 4 .  • PH  !  •>; 

EQWRITE ; 

END 

ELSE 

BECIN 

WRITE (  FF  ,  *  '.MliS,*  !  E  =  '  , EN : 8 ; 4 ) ; 
WRITE  <  FF  ,  ‘  '  ,  HE  ;  7  :  4  ,  '  PH  !  '); 

EQWRITE ; 

END; 

END; 

FOR  KK ; * 1  TO  80  DO  MATC  KK3  ; * '  •; 

FOR  KK:-1  TO  15  DO  NAME C 1 , KK 1 ; -TEMP C KK 3 ; 

END; 

PROCEDURE  EQB4; 

BEGIN 

FOR  KK-l  TO  80  DO  MATCKK2 : - '  •; 

1 1  -1  ; 

FOR  J:«l  TO  I  DO 
BEGIN 

IF  JOl  THEN 
BEGIN 

MATC III; -SC  J3 ; 

I  I : -I  I ♦  l  ; 

MATCH!  ; 

I I  :  -I  I  4-1  ; 

END; 

IF  CCJIOI  THEN 
BEGIN 

GETVAL (  CC  J 1 , SAVE)  ; 

FOR  KK : ■ 1  TO  15  DO 

IF  SAVE  C  KK 1 <  > '  •  THEN 

BEGIN 

MATC III: -SAVECKK3 ; 

II  : ■ I  I  4  1  ; 

END; 

MATC I  I  3  : - ’ *  ’  ; 

I  I  :  -  I  I  4  l  ; 

END  ; 

FOR  KK-l  TO  13  DO 
SAVE  C  KK 1  : -XC J , KK3 ; 

FOR  KK : - 1  TO  15  DO 

IF  SAVE  C  KK  3  <  >  '  1  THEN 

BEGIN 

MATC III; -SAVE  C  KK  3 ; 

II : - 1 1  4 1  ; 

END; 

MATC III:-1  • ; 

II  : -  I  I  4  1  ; 

I  I  ;  -  I  I  4 1  ; 

END, 

MAT! I  I  1  :  - ' • '  ; 

END; 


PROCEDURE  LOCSUM ; 

BEGIN 
I  :  >0  , 

SUMAIO  =0  0 ; 

FOR  KK  =1  TO  NION  DO 
BECIN 

SUMAIO:  -SUMAIO+COEFCKK+23*  LOG ( A  I ONC  KK I > ; 
END  , 

END; 


PROCEDURE  EQB 1  , 

BEGIN 

MI  : -Ml +1  ; 

IONL INCMI 3 : -0 ; 

I  F  <  AM  CM  3  <  >  1  )  OR  (AMCAIOl)  THEN  I ONL INC M I  3  : - 2 ; 

I  F  (  AMCM3  <  >  1  >  AND  (AMCAIOl)  THEN  I ONL  I NC  M I  1  :  -  1  ; 

IF  MI-INDEX  THEN 
BEC  IN 

IF  INDEX  <  > 1  THEN  DRAWERC  79 )  ; 

INDEX  :  -  I NDE X 4-7  0  ; 

LNFD  < 15); 

WRITELH(CHR< 1 2  > )  ; 

DRAWER ( 79  >  ; 

VRITELNC FF ,  ’  !  !  '  >  ; 

WRITE(FF, 'NUMBER  !  EQUILIBRIUM  !'); 

WR ITELN ( FF  ,  '  BALANCED  CHEMICAL  EQUATION'); 

WRITELNt  FF ,  '  •  !  '  )  ; 

DRAWER ( 7  9 >  ; 

END, 

ACT 1  -AMCMI ; 

ACT2 : -AMC AI ; 

FOR  HH.-l  TO  15  DO  NAME  CNI ON +  3 , HH3  : -SP  C  A , HH 1 ; 

COEFCNI ON* 3  I  : -METALCM1 ; 

COEFC 1 1 . -METAL CA1 ; 

FOR  I ; -1  TO  NION  DO 

COEFC 1  +  2  I : -COE F CNI ON* 3  I  * IONC I , A J -COE F C 13 « I ONC I  , Ml  ; 
SUMNUM ; 

NUMB : -NION+5 ; 

FOR  HH : - 1  TO  15  DO  TEMP C HH 1 ; -NAME C 1 , HH 1 ; 

LINEC 1  .MI  1  : -M; 

LINEC  2 .MI  1  : -A; 

IF  ( <  COE F  CNI ON*  51-0) AND ( COEFCNI ON*  4  I <  0 ) ) OR ( COEF CNI ON+5 3 <0 
UPDATE ; 

END; 


PROCEDURE  EQBAL ANCER 
BEGIN 
MI  ;  -0  ; 

NAME  C  2 .  1 1  : - ' H ' 
NAME C 2 . 23 ; - ' 2 ' 
NAME  C  2 , 33 : -'O' 
NAME  CNION+4 . 1 3 
NAME  CNI ON+  4.23 
NAMECNION+S , 1 3 
NAME CNI ON +5 , 2  3 


NPP : -NP-1 ; 

FOR  M: -1  TO  NPP  DO 
BECIN 

FOR  HH : ■ 1  TO  IS  DO  NAMEC 1 , HH3 ; 
K ; -M+l ; 

FOR  A: -K  TO  NP  DO 
BECIN 


EQB  1  ; 
LOCSUM; 
FINDSIGN; 
EQB 4  ; 
CALCEQEQ ; 
END; 


END; 


-SPCM.KH3 ; 


)  THEN 


END; 


PROCEDURE  DIACRAM, 

BEGIN 

DRAWER  <78 )  ; 

IONS : >0 ; 

NON I  ON :  =0 ; 

FOR  I : ml  TO  NP  DO 

IF  AMCIJ-1  THEN  NON I  ON :  =NONION+ 1 
ELSE  IONS : - 1 ONS* 1 ; 

IF  NONION-O  THEN 
BEGIN 

LNFDC 10); 

WRIT£LN<  FF ,  '  **  NO  SOLID  SPECIES  DIAGRAM  **•); 

END 

ELSE 

BECIN 

IONDIG : =0 ; 

STABLE  1  ; 

END  ; 

IF  I ONS  <  » 1  THEN 
BECIN 

LNFD< 10); 

WRITELN< FF , '  «**  NO  DISSOLVED  SPECIES  DIAGRAM  •««•> 

END 
ELSE 
BECIN 

I OND I G : = 1 ; 

STABLE  1  ; 

END; 

END; 


■T'  )  AND  ( OPT 1  t  >  1  P 1  >  DO 


PROCEDURE  OUTMENU; 

BEGIN 

WHILE  lOPTIO 
BEGIN 

WRITELN( CHR< 12) ) 

IF  AB- FALSE 

THEN  WR ITELN< ' 

ELSE  WRITELN< 1 
LF  <  3 )  ; WRITELNt  • 

LF  < 1 )  ;WRITELN( * 

LF  <  2  >  ; WRITE  <  1 
RE ADLN< OPT1 > ; 

WRITELNI CHR< 1 2  > )  ; WRITE LN; 

END; 

IF  OPT 1 «  1  P  1  THEN  FILENAME ; -'PRINTER 
ELSE  FILENAME CONSOLE ; 

END; 


DATA  OUTPUT  OPTIONS  1 ) 
EQUATION  LISTING  OPTIONS1 > 


P  -  PRINTER1 > ; 
T  -  TERMINAL  1  ) 
ENTER  OPTION: 1 > 


( *  *N+  «  ) 

BEGIN 
R  E  P  E  AT 

WRITELN! CHR! 12)); 

INITIALIZE ; 

DATAINPUT; 

METAL  INPUT ; 

WRITE! '  MODIFY  DATA  »  CY/NJ:  •>; 

RE ADLN( OPT1 > ; 

IF  OPT 1  -  '  Y  '  THEN  DATACORRECT; 

WRITELN! CHR( 12) > ; 

RE  PEAT 

IF  NO PT  =  ‘  Y  '  THEN  DATACORRECT; 

RE  P  E  AT 

LF ( 1 > ; WRITE ( '  DISPLAY  TABULATED  DATA  »  CY/NJ 
READLNtOPTl ) ; 

AB : -FALSE ; 

IF  OPTl-'Y'  THEN 
BECIN 

OUTMENU ; 

REWRITE  <  FFF .FILENAME  >  ; 

NONMETALDI SPLAY ; 

METAL D I  SPLAY ; 

CLOSE! FFF ) ; 

OPT l  : -  *  1  ; 

WHILE  (  OPT  1  <  >  '  Y  *  )  AND  !OPTl!)'N'>  DO 
BEGIN 

WRITE!'  MODIFY  DATA  »  CY/NJ:  '); 
REAOLN! OPT1 ) ; 

END; 

IF  OPTl-'Y'  THEN  DATACORRECT; 

END; 

AB : -TRUE ; 

OUTMENU ; 

REWRITE !FF, FILENAME) ; 

EQBALANCER ; 

DIAGRAM; 

CLOSE ! FF ) ; 

LF  !  2  )  ; 

WRITELN! •  EXECUTION  COMPLETE ....'); 

LF ! 3  )  ; 

WRITE!'  EXECUTE  AGAIN  ?  CY/NJ:  '); 

READLN! NOPT) ; 

IF  NOPT-'N'  THEN  EX  IT! PROGRAM)  ; 

INDEX : -1 ; 

WR ITELN! CHR ! 12) ) ; 

LF!  2  >  ; 

WRITE! '  REPEAT  WITH  SAME  DATA  ?  CY/NJ:  '); 
READLN! NOPT) ; 

UNTIL  NOPT-'N’ ; 

LF ( 2  >  ; 

WRITE! '  MODIFY  CURRENT  DATA  SET  »  CY/NJ:  '); 
READLN! NOPT) ; 

UNTIL  NOPT-'N' ; 

UNTIL  NOPT- ' X ' ; 

END 


(  *  <s+  *  > 


UNIT  CETVALUE.  INTRINSIC  CODE  23  DATA  24; 

INTERFACE 

TYPE 

AO  t  5* ARRAY  Cl  1 S  3  OF  CHAR; 
AOI15.ARRAYC1  15!  OF  INTEGER; 

PROCEDURE  GETVALC X : INTEGER ; 

VAR  CCHAR : A013 ) ; 

IMPLEMENTATION 

VAR  HH , I  . INTEGER ; 

I  NT  : AO  I  1 5 ; 

VAL  : AO  IIS; 

PROCEDURE  GETVAL ; 

BEGIN 

FOR  I ; - l  TO  15  DO 
BEGIN 

INTC I  3  :  =  0  ; 

VALC I  3  ; -0  ; 

CCHARC 13:='  * ; 

END; 

HH;  -1  ; 

WHILE  X>«  1  DO 
BEGIN 

INTCHH3 : -X ; 

X : .X  DIV  10 ; 

HH: «HH+l ; 

END; 

VALC  HH  3 : ■ INTCHH3  ; 

FOR  I : -1  TO  (HH-l )  DO 

VAL  C  HH+ 1-13 : -INTC 13-10* INTC 1  +  13 
FOR  I.-l  TO  HH  DO 

CCHARC I  3 : «CHRC VALC I 3+ORDC  ’O')); 
FOR  I: -l  TO  <  HH- 1 )  DO 
CCHARC I  3 : -CCHARC 1  +  13; 

CCHARCHH3  ; -  *  '  ; 

END  ; 

BECIN 

END 


C  *  SS*  •  ) 


UNIT  CONVERT.  INTRINSIC  CODE  16  DATA  17; 

INTERFACE 

TYPE 

AOCH3 * ARRAY C 1  .31  OF  CHAR; 

AOCH1 t-ARRAYCl .  11]  OF  CHAR; 

VAR 

NUMBER  INTECER; 

CH3  AOCH3; 

ERROR  BOOLEAN; 

CH11  AOCHll; 

RNUMBER : REAL ; 

PROCEDURE  CHARTOINT; 

PROCEDURE  CHARTOREAL, 

IMPLEMENTATION 

VAR  I  INTECER; 

BLANX : BOOLEAN; 

POWER : REAL ; 

PROCEDURE  CHARTOINT; 

BECIN 
I  :  *  3  ; 

NUMBER : >0 ; 

BLANK-FALSE ; 

POWER: -1 . 0 ; 

ERROR: -FALSE; 

WHILE  <  CH3  C  I  J  -  '  ')  AND  <  X  >  1  >  DO 

IF  <  I  -  1  >  AND  ( CK3  C I  1 « ‘  •  > 

THEN 

BECIN 

WR I TE  LN ' 

WRITELNt 'NUMBER  EXPECTED*); 

ERROR : -TRUE ; 

END 

ELSE 

BECIN 

WHILE  (ERROR-FALSE)  AND  (I>0)  AND  (BLANK-FALSE) 
IF  CH3CIJ  IN  I  ’ 0 •  .  .  * ? *  .  *  -  * ]  THEN 
BECIN 

IF  CH3  C I J - ' -  *  THEN 

NUMBER: --NUMBER 

ELSE  NUMBER : -NUMBER* ( 0RD( CH3  C  l  J  ) -ORD(  *0*  > 
POWER: -POWER* 10 . 0 ; 

I  :  -  I  -  1  ; 

END 

ELSE 

IF  CK3CI1-*  *  THEN 
BLANK : -TRUE 
ELSE 
BECIN 

WRITELN; 

WR ITELN ( '  IMPROPER  NUMBER' >; 

ERROR : -TRUE ; 

END; 


DO 


*TRUNC( POWER) 


END; 


END; 


PROCEDURE  CHARTOREAL. 

BEGIN 
I  -  1  1  . 

RNUMBER  *  0 , 

BLANK  .FALSE ; 

POWER  =10; 

ERROR  .FALSE ; 

WHILE  ( CH l l [ I]« '  •)  AND  <  I  >  1  >  DO  I  =1-1; 

IF  (1  =  1)  AND  <  CH 11(1}='  •>  THEN 

BECIN 

WRITELN ; 

WRITELNI*  NUMBER  EXPECTED'); 

ERROR  -TRUE ; 

END 

ELSE 

BEGIN 

WHILE  <  ERROR. FALSE  >  AND  <  I  >  0  >  AND  (BLANK-FALSE)  DO 
IF  CHlttll  IN  C  ‘ 0  *  .  .  ' 9 '  ,  *  - '  ,  '  .  •  3  THEN 
IF  CHI  1 t I  I  =  '  .  '  THEN 
BEGIN 

RNUMBER  -RNUMBER /POWER ; 

POWER : -V . 0 ; 

I  : -I-  1  ; 

END 

ELSE 

BEGIN 

IF  CHllCl!-'-'  THEN 
RNUMBER : --RNUMBER 
ELSE 
BEGIN 

RNUMBER ; -RNUMBER* ( ORD < CHI l t I  3 ) -ORD<  ' 0 ‘  ) > ‘POWER; 
POWER; -POWER* 10 . 0 ; 

END; 

I  :  *  I  -  1  ; 

END 

ELSE 

IF  CHI  It  II-'  '  THEN 
BLANK; -TRUE 
ELSE 
BECIN 

WRITELN ; 

WRITELN< '  IMPROPER  NUMBER'  >  ; 

ERROR; -TRUE; 

END; 


END; 

BECIN 

END 


END; 


UNIT  STABLE,  INTRINSIC  CODE  26  DATA  27; 


INTERFACE 

USES  CONVERT, INDATA , 

TYPE 

AOINT2X395-ARRAYC 1  2,1  595 ] OF  INTECER; 

AO I NTS  9  5  ■ ARRAY  C  1  595  1  OF  INTECER, 

AOR2X595  -ARRAYC1  2,1  59S30F  REAL; 

VAR  FF : TEXT ; 


F I LENAME : STRINC ; 

MI  : INTECER ; 

LINE  : AO  I NT2  X  5  9  5  ; 
LTYPE  : AO  I NTS  9  5 ; 
EQUAT  AOR2X595; 
IONDIG  ; INTECER; 
IONLIN  : AO  I NTS  9  5 ; 


PROCEDURE  STABLE  1  ; 

PROCEDURE  LNFDIQQ : INTECER) ; 
PROCEDURE  DRAWERIQQ : INTEGER ) ; 


IMPLEMENTATION 

VAR 

NUMB  1  , LINNO, LINNN, LTYPE 1 

II ,L1 ,L2, INDEX , J INDEX , I , N,NN,PASS 

PHMIND , PHMAXD , K , J J . HH , J 

LINE1 

EQUAT 1 

LOVER, UPPER 

PHTEST 

INTERC , EMINU , EMINL , PHI , PH2 
ABC, POT1 . POT2 , POT 1 1 , POT22 , El . E2 
PHMIN, PHMAX, EMAXU , EMAXL , PH I NT 
FLAG  1  , FLAG 3 . FL AG4 . FLAG 6 , FLAGGER 
FLAGS , FLAC7 , FLAGS , FLAG9 , FLAG  10 


ARRAY  Cl  .  .34]  OF  INTEGER; 
INTECER; 

INTECER; 

ARRAY  Cl.. 2,1.. 341  OF  INTEGER 
ARRAY  Cl  .  .2,1.  .34]  OF  REAL; 
ARRAYC1 .. 4 , 1 . . 343  OF  REAL; 
ARRAY  C 1  .  . 35 1  OF  REAL ; 

REAL  ; 

REAL  ; 

REAL; 

BOOLEAN; 

BOOLEAN; 


PROCEDURE  LNFD; 

VAR  Q; INTEGER; 

BECIN 

FOR  0;-l  TO  QQ  DO  VRITELN! FF ) ; 
END; 


PROCEDURE  DRAWER; 
VAR  Q  INTEGER; 
BEGIN 


FOR  Q;-l  TO  QQ  DO  WRITE! FF 
WRITELN(FF) ; 

END; 


PROCEDURE  RITER( ARC  1  , ARC 2 ; REAL ; 

ARG3 , ARG4 : REAL ; 

ARCS ; INTEGER) ; 

BEGIN 

WRITE (FF.*  I  * , ARC  1:7:3,'  !  ,,ARG2;6:2); 

WRITE  <  FF ,  '  I  '  , ARG3 : 7 : 3 , 1  I  '); 

WR I TE  LN  <  F  F , ARC  4:6:2,'  !  '  , ARCS : 3 ,  '  !'); 

END; 


PROCEDURE  SUB 1 ; 

BEGIN 

IF  <  POT1  1  <  POT1  >  AND  <PHMINDOO>  THEN 

R ITER  <  POT1 , UPPERI 3 , 1 ] , POTl 1 , LOVERC3 , 1 ] , NUMB 1 C PHMIND 1 > ; 
IF  ( POTl 1< POTl)  AND  (PKMIND-0)  THEN 
BECIN 

WR ITE ( FF , '  !  * • ,POTI :7: 3, '  !  '); 

WRITE !  FF  , UPPERC  3,11:6:2,'  !  *  , POTl l : 7 : 3 ) ; 

WR ITELN( FF , ‘  !  • , LOVERC 3 , 1 3 . 6 : 2 , '  !  !'); 

END; 

IF  <  POT2  2  <  POT2 )  AND  (PHMAXD! >0)  THEN 

R ITER (POT2 .UPPERC 4 ,N1 , POT2  2 , LOWER C  4 ,NN1 , NUMB 1 C PHMAXD 1 > ; 
IF  ( POT2  2 ! POT2 )  AND  (PHMAXD-0)  THEN 
BEGIN 

WRITE ! FF ,  '  1  '  , POT2 : 7 : 3 ,  '  !  •>; 

WRITE ( FF , UPPERC  4 , N1 ; i ; 2 , '  !  '  , POT22 : 7 : 3 ) ; 

VRITELN! FF , '  \  ‘ , LOVERC 4 ,NN J : 6 : 2 , '  I  !'); 

END; 


END; 


PROCEDURE  STABLE  1 1  , 

BECIN 

FOR  I  :  -1  TO  N  DO 
BECIN 

POT1  :  -UPPER t  1,11  ♦UPPERC  2  ,  I  3 "UPPERC  3,13; 

POT2 : -UPPERC 1 , I  3 ♦UPPER t  2 , I  3 ‘UPPER C  4  ,13; 

R I TER <  POTl  , UPPERC  3 ,  I  3  , POT2 , UPPER C 4,13 , LINNOC 13); 

END; 

FOR  I : -1  TO  NN  DO 
BEGIN 

POTl : -LOWER! 1  , I  3 ♦LOVERC  2 ,  I  3 "LOVERC 3. 13 ; 

POT2 ; -LOVERC 1.13 ♦LOWERC2 , I  3 "LOVERC 4 ,13; 

RITER(  POTl , LOVERC  3,13, POT2 , LOVERC  4 ,13, LINNNC 13); 

END; 

POTl  . -UPPERC 1 , t  3 ♦UPPERC 2 . 1 3 "UPPERC 3 ,13; 

POTl 1  r -LOVERC 1 ,  l 3  +  LOVERC 2 . 1 3  *LOVERC  3 ,13; 

POTl 1 : - POT 11+0 . 0001 ; 

POT2 : -UPPERC 1 ,N3 ♦UPPERC 2 ,N3 "UPPERC 4 ,N3 ; 

POT2  2 ; -LOVERC 1 ,NN3 ♦LOVERC 2 , NN3  * LOVER C  4 , NN3 ; 

POT22 : -POT22+0 . 0001 ; 

SUB1  ; 

FLACCER: -TRUE ; 

END; 

PROCEDURE  STABLE 1 0 ; 

BECIN 

LNFD ( 2 ) ; 

VRITE  <  PF ,  *  *  ) ; 

DRAVER <  5  7)  * 

VRITE (  FF  ,  *  '  1  POTENTIAL  11  PH  1  !'); 

VRITELNC  FF ,  *  POTENTIAL  2  1  PH  2  I  LINE  NO.  I'); 

VRITE  <  FF ,  *  *); 

DRAVER (37) ; 

END; 

PROCEDURE  STABLE? ; 

BECIN 

FLAG 1 0 : -FALSE ; 

IF  EMAXL  > -EMAXU  THEN 
BECIN 

VHILE  FLAG 1 0-FALSE  DO 
BECIN 

IF  (UPPERC4.N1>  PHTESTC  K3 ) AND ( N> 1 >  THEN  N : -N- 1 ; 

IF  (LOVERC4 , NN3 > PHTESTC K 3 ) AND(NN> 1 )  THEN  NN ; -NN- l ; 

K  '  — K—  1  ' 

EMAXL; -LOVERC 1 .NN3 ♦LOVERC 2 . NN 1* PHTESTC K 3 ; 

EMAXU: -UPPERC 1 ,N3^UPPERC2 ,N3 "PHTESTCK3 ; 

IF  EMAXU) -EMAXL  THEN  FLAC1 0 : -TRUE ; 

END; 

PHINT: -(UPPERC 1 ,N3 -LOVERC 1 ,NN3  > / ( LOVERC 2 ,NN3 -UPPERC 2 ,N3 > 
LOVERC 4 ,NN3 : -PHINT; 

UPPERC4 , N3 ; -PHINT; 

END; 

END; 

PROCEDURE  STABLES ; 

BECIN 

IF  FLACS-FALSE  THEN 
BECIN 

PHINT; -(UPPERC 1 , 1 3 -LOVERC 1 , 1 1 >  M LOVERC 2 , 1 3 -UPPERC 2 , 1 3 ) ; 
LOVERC  3 , 1 3 ; -PHINT; 

UPPERC  3.13 ; -PHINT; 

END; 

END; 

PROCEDURE  SUB4CL: INTEGER) ; 

BECIN 

LOVERC 1 , L  3 : -LOVERC 1 .L^ll ; 

LOVERC 2 ,L3 : -LOVERC 2 . L+ 1 3 ; 

LOVERC 3, LI : -LOVERC 3, L+ 13 ; 

LOVERC  4 . L3 : -LOVERC 4 ,L^ll ; 

LINNNC L 3  ;  -LINNNC L*1  3  ; 

END; 

PROCEDURE  SUBS (L: INTEGER) ; 

BEGIN 

UPPERC  l  .LI  :  -UPPERC  1  ,L*11  ; 

UPPERC 2. L 3 : -UPPERC Z.L* 11 ; 

UPPERC 3, L3 .-UPPERC 3, L^l 3; 

UPPERC 4. L 3 : -UPPERC 4 ,L*11 ; 

LINNOC  LI :  -LINNOC  LM  3  ; 

END; 


PROCEDURE  STABLE?; 

VAR  L : INTEGER; 

BEGIN 

FLAGS ; -FALSE ; 

IF  EMINL  <  EMINU  THEN  FL AG8 : -TRUE ; 

IF  FLAC8-FALSE  THEN 
BECIN 
I  :  -  l  ; 

FLAG? : -FALSE ; 

WHILE  FLAC9 -FALSE  DO 
BECIN 
I  : -  I +t  ; 

ABC  *  — PHTESTC I  3  ' 

IF  (ABC<0 . 00000000  1  )  AND  < ABC > -0 . 0 0 0 0 0 0 0 1  >  THEN  ABC:-0.0; 
EMINL: -LOVERC 1 ,  1 1 +LOWERC  2 . 1 3* ABC; 

EMINU: -UPPER 1 1 . I  3  +UPPERC 2 , 13*ABC; 

I F ( EM I NU < - EM I NL  >  AND  (I-K)  THEN  FLAC7 : -TRUE ; 

IF  ( FLAG7 -FALSE )  THEN 

IF  EMINL  <  EMINU  THEN  FLAG9 : -TRUE ; 

IF  <  FLAG? -FALSE  >  AND  ( FLAG9-FALSE  >  THEN 
BECIN 

IF  LOWERC  4  ,  1  3  < -PHTESTC I  3  THEN  LINNNC13:«0; 

IF  UPPERC4 , 1 3 <«PHTESTC I  3  THEN  L INNOt 1 3 : -0 ; 

IF  LINNNC 13-0  THEN 
BECIN 

NN : -NN- 1 ; 

FOR  L : - 1  TO  NN  DO  SUB4(L>; 

END; 

IF  L I NNO  C 1 3  - 0  THEN 
BEGIN 
N  :  -N- 1  ; 

FOR  L : - 1  TO  N  DO  SUBS(L); 

END; 

END; 

IF  I -35  THEN  FLAG9 : -TRUE ; 

END; 

END; 

END; 

PROCEDURE  SUB3 ; 

BECIN 

PHTESTC K1 : -UPPERC4  .  I  3  ; 

IF  < I -N)  AND  (JJ-NN)  THEN  FLAGS : -TRUE ; 

IF  FLAGS-FALSE  THEN 
BEGIN 

IF  I <N  THEN  I :  -1  +  1 ; 

IF  JJ <NN  THEN  JJ.-JJ+l; 

K : — K+ 1 ; 

END; 

END; 


PROCEDURE  STABLES ; 

BECIN 

IF  FLAGS-FALSE  THEN 
BEGIN 
I  :-l  J 
JJ : -I  ; 

K:-2; 

PHTESTC 1 3 : -PHMIN; 

FLAGS  -.-FALSE; 

WHILE  FLAGS-FALSE  DO 
BEGIN 

IF  UPPERC4, 13 (LOWERC 4, JJ3  THEN 
BEGIN 

PHTESTC K 3 :-UPPERC4, 13  ; 

K : -K+ 1 ; 

I : -I  +  l ; 

END 

ELSE 

IF  LOWERC4, JJ3<UPPERC4, 13  THEN 
BECIN 

PKTESTCK3 : -LOWERC 4 , J J 3 ; 

K : -K+ 1 ; 

JJ : -JJ*  1  ; 

END 

ELSE 

IF  UPPER C  4,1 3 -LOVERC  4 , JJ3  THEN  SUBS; 

END; 

END; 

END; 


i 


PROCEDURE  STABLES; 

BEGIN 

EMINL : - LOWER £  1 , 1 3  4- LOWERE  2 ,  1 3 ‘PHMIN; 

EMINU: -UPPER C 1 .13  4-UPPERE2.1 3 ‘PHMIN; 

EMAXL  ;  ■  LOWER  £  1  ,NN3  4 LOWE RE  2  ,  NN  3  *  PHMAX  ; 

EMAXU  :  -UPPER  £  l  ,NUUPPERt2  ,N3*  PHMAX; 

IP  <  EMINL <  EM INU )  AND  < EMAX L <> EMAXU >  THEN  FLAGS : -TRUE ; 
END  ; 

PROCEDURE  SUBS; 

BEGIN 

LNFD ( 2 > ; 

WRITE ( FF  ,  '  ’  >  ; 

nniup p / s  s  \  ■ 

WRITE  <  FF  ,  ‘  !  POTENTIAL  l  !  PH  1  !’>; 

WRITELN( FF ,  '  POTENTIAL  2  1  PH  2  !  LINE  NO.  I  *  >  ; 
WRITE<  FF .  •  ’  >  ; 

DRAWER ( 35 ) ; 

IF  NOO  THEN 

FOR  I : -1  TO  N  DO 
BEGIN 

POT1 : -UPPER! 1 . I  3  +UPPERE 2 ,  I3*UPPERE3,  II  ; 

POT2 : -UPPER E 1  ,  I  1 +UPPERE  2 ,  I 3‘UPPERE  4 ,13; 

R ITER ( POT 1 .UPPER! 3, 13 . POT2 , UP PER E 4 , 13 . LINNOEI3) ; 
END 

ELSE 

FOR  l : -t  TO  NN  DO 
BEGIN 

POT1  :  -LOWER!  1  ,  I  3  4-LOWERE  2  ,  I  3  *  LOWE  RE  3  ,13; 

POT2 : -LOWERE 1 ,  I  3  4 LOWER! 2 ,  I  3  *  LOWER  E  4 ,  I  3  ; 
RITEREPOTl  , LOWERE 3 , I  3  , POT2 , LOWERE  4.13 , LINNNE 13); 
END; 

END; 


PROCEDURE  INTER IMWRITER; 

BEGIN 

IF  (N-0)  OR  (NN-0 )  THEN 
BEGIN 

IF  N-0  THEN 
BEGIN 

LNFDEZ) ; 

WRITELNE FF , *  NO  UPPER  LIMIT’ >; 

END; 

IF  NN- 0  THEN 
BEGIN 

LNFDE2); 

WRITELNEFF.’  NO  LOWER  LIMIT’); 

END; 

IF  ENOO)  OR  ENNE  >  0  )  THEN 
SUBS  ; 

IF  PHMIN) -2.0  THEN 
BEGIN 

LNFDE I > ; 

WRITE E  F F  ’  '  )  • 

WRITEE FF 1 ’MINIMUM  PH  -  ’ , PHMIN : 6 : 2 , ’  LINE  NO.  •>; 
WRITEINE FF ,NUMB1 EPHMIND1 :4)  ; 

END; 

IF  PHMIN! --2 . 0  THEN 

WRITELNEFF,'  MINIMUM  PH  -  ’  . PHMIN : 4 ; 2 ) ; 

IF  PHMAX < IS . 0  THEN 
BEGIN 

WRITEEFF,'  MAXIMUM  PH  -  ’ , PHMAX : 4 : 2 ) ; 

WRITELNE ’ PHMAXD-  ’,PHMAXD:4); 

WRITELNEFF,'  LINE  NO.  ’  , NUMB I E  PHMAXD 3 : 4 ) ; 

END; 

IF  PHMAX)- 14.0  THEN 

WRITELNEFF,'  MAXIMUM  PH  -  ' , PHMAX : 4 : 2 ) ; 

FLAGGER; -TRUE; 

END; 

END; 


PROCEDURE  TEMP2 ; 

BEGIN 

FLAG  3 : -FALSE ; 

NN : -NN  + 1 ; 

LOVERC 1  , NN  3  : -EOUAT1 C 1 . LI  3 ; 

LOVERC  2 , NN I  : -EQUAT1 C  2 , L 1  I  ; 

LOVERC  3 ,NN3  : -PHI  ; 

L INNNCNN3  : -NUMB  1 CL1  3  ; 

LOVERC  4 ,  NN3 . -PHMAX ; 

IF  PASS < 1  THEN  FLAC3 : -TRUE ; 

IF  FLAG3-F ALSE  THEN 
BEGIN 

LOVERC  4 ,NN3  :  -PH2  ; 

IF  PH2  > - PHMAX  THEN  LOVER C 4 , NN 3  : - PHMAX ; 
IF  PH2> -PHMAX  THEN  FLAG3 : -TRUE ; 

IF  FLAG3-FALSE  THEN 
BEGIN 

PHI . -PH2 ; 

El  : -EQUAT1 C 1  . L 2  3  + EQUATl C 2 , L2  3* PHI ; 
LI  :  -L2 ; 

END  ; 

END; 

END; 

PROCEDURE  TEMPI ; 

BEGIN 

FLAG  1 : -FALSE ; 

N : -N+ 1 ; 

UPPERC l ,N3 ; -EQUAT1 C 1 ,L1 3 ; 

UPPERC2.N3  : -EQUAT1 C  2 , L 1 3 ; 

UPPERC  3 , N3 ; -PHI ; 

L INNOCN3 : -NUMB  1 C L 1 3 ; 

UPPERC  4 , N3 : -PHMAX ; 

IF  PASS < 1  THEN  FLAG 1  : -TRUE ; 

IF  FLAG l -F ALSE  THEN 
BEGIN 

UPPERC  4 , N3 : -PH2 ; 

IF  PHI > -PHMAX  THEN  UPPERC 4 ,N3 : -PHMAX ; 
IF  PH2  > -PHMAX  THEN  FLAG l ; -TRUE ; 

IF  FLAC 1 -FALSE  THEN 
BEGIN 

PH  t  *  * PH2 * 

El  : -EQUATl C 1  , L2  3 ♦ EQUATl C2,L23*PH1  ; 
LI : -L2 ; 

END; 

END; 

END; 

PROCEDURE  LOVBNDLNS ; 

BEGIN 

FOR  I : -1  TO  II  DO 

I F ( LINE 1C2|I3-J)  AND  <  LTYPE 1 C I  3  - 2 )  THEN 
BECIN 

PASS  :  -  PASS  4- 1  ; 

IF  PASS-1  THEN 
BECIN 

El  .-EQUATl  Cl  .  U4-EQUAT1C2  ,  I  3  *  PHI  ; 
LI  :  -I ; 

END 

ELSE 

BEGIN 

E2 : -EQUATl C 1 , I  3  +EQUAT1 C  2 , I  3  *PH1 ; 
IF  E2  > -E 1  THEN 
BEGIN 
El : -E2  ; 

L  1  :  -  I  ; 

END; 

END; 

END; 

END; 


PROCEDURE  I NTPTSLO ; 

BECIN 

PASS : -0; 

FOR  I  : - 1  TO  II  DO 
BECIN 

IF  (LINE1 C 2 , I 3-J >  AND  < LTY P E 1 C I  3  -  2 >  AND 
( EQUAT1 1  2 ,  I ) >  EQUAT1 C  2 , L 1  I >  THEN 
BECIN 

PASS : -PASS ♦ 1  ; 

IF  PASS  <-  1  THEN 

BEGIN 

PH2 :  -EQUAT1 I  1  ,  LI 3-EQUATl C l . 13 ; 

PH2 : -PH2 / <  EQUAT1 C  2 , I 3-EQUAT1 C2 . LI  3  >  ; 

L  2  :  -  I  ; 

END 

ELSE 

BECIN 

INTERC : -EQUATt  C 1 , LI 3-EQUAT1 C 1 , I  3 ; 

INTERC : - INTERC / ( EQUAT1 C2 ,  I 3-EQUATl C2 . LI  3  > 
IF  < INTERC<-PH2 >  AND  ( INTERO -PHI )  THEN 
BECIN 
L2 : -I  ; 

PH2 : = INTERC 
END  ; 

END; 

END; 

END; 

TEMP2 ; 

END; 

PROCEDURE  INTPTSH I ; 

BEGIN 

PASS ; -0  ; 

FOR  I:  -1  TO  II  DO 

IF  (LINE1C1 . I3-J)  AND  ( LTYPE 1 C  I  3  -  2  )  THEN 
IF  EQU  AT 1  C  2  ,  I  3  <  EQU AT 1  I  2 , L 1 3  THEN 
BECIN 

PASS: -PASS+1 ; 

IF  PASS  <  «  1  THEN 
BECIN 

PH2 : -EQU ATI C 1 , L 1 3 -EQUAT1 C 1 ,13; 

PH2 : -PH2 / <  EQUAT1 C  2 , I  3 -EQUAT1 C  2 , L 1 3 )  ; 

L2 : -I  ; 

END 

ELSE 

BECIN 

INTERC : -EQU ATI [ 1 , LI 3-EQUATl C 1 , I  3 ; 

INTERC : m INTERC / ( EQU ATI C 2 , I  3 -EQUAT1 C 2 , L 1 3 ) ; 
IF  INTERC <  «PH2  THEN 
BECIN 
L2 : -I  ; 

PH2 : ■ INTERC ; 

END; 

END; 

END; 

TEMP1 ; 

END; 

PROCEDURE  UPBNDLNS ; 

BECIN 

PASS . -0  ; 

FOR  I : -1  TO  II  DO 

IFILINEl C 1  , I 1-J)  AND  < LTYPE 1 C I  3 -2 >  THEN 
BECIN 

PASS : «  P ASS  +  1  ; 

IF  PASS-1  THEN 
BECIN 

El : -EQU ATI C 1 ,  I  3  +EQUAT1 C  2 , I3*PH1 ; 

LI  :  -I ; 

END 

ELSE 

BECIN 

E2 : -EQU ATI C 1 , I  3  + EQU ATI  I  2 , I  3  *PHl ; 

IF  E2  < -E 1  THEN 
BECIN 
El : -E2 ; 

LI  :  -I  ; 

END; 

END; 

END; 

END; 


PROCEDURE  STABLE*; 

BECIN 

PHMAX : =1 6 . 0 ; 

PHMAXD : ■ 0 ; 

FOR  I : *1  TO  II  DO 

IF  LINE  1  I  1  .  I  3 -J  THEN 
IF  LTYPE1CI3-1  THEN 
IF  EQUAT1 C  1  .  I  3  <  PHMAX  THEN 

“SJSuC.O.AT.,,...; 

PHMAXD  ; -  I ; 

END; 

END; 

PROCEDURE  STABLE3 ; 

BEGIN 

PHMIN; --2 . 0 ; 

PHMIND : -0 ; 

FOR  I : -1  TO  II  DO 

IF  ( L INE 1  C  2  ,  I  3  ■  J  >  THEN 
IF  ( LTYPE 1 113-1)  THEN 

IF  (EQUAT1 C 1 ,  I  3 >PHMIN>  THEN 
BCG  IN 

PHMIN: -EOUAT1 C  t ,  I  3 ; 
PHMIND: -I ; 

END  ; 

END  ; 

PROCEDURE  STABLE2; 

BEGIN 

II  :»I  14-1  ; 

L INE 1 C 1,113; -LINEC 1 ,13; 

LINE1 12 , 113: -LINEC2 . I  3 ; 

LTYPE 1 C 113: -LTYPEC I  3  ; 

EQUATl C 1,113; -EOUATC 1,13; 

EQUAT1 1  2 , 113: -EOUATC 2 , I  3 ; 

NUMB  1 C I  1 3 : *  I ; 

END; 

PROCEDURE  START; 

BEGIN 

INDEX  ;  -1  ; 

JINDEX : -0  ; 

END; 

PROCEDURE  SETFLAGS; 

BECIN 

FLAGCER: -FALSE  ; 

FLAGl : -FALSE ;  FLAG3 : -FALSE ; 

FLAC4 : -FALSE ;  F L AG 5 : - F AL SE ; 

FLAG* : -FALSE ;  FLAG7 ; -FALSE ; 

FLAGS ; -FALSE ;  FLAG? : -FALSE ; 

FLAG  1 0 : -FALSE; 

END; 


PROCEDURE  STABLE1 ; 

BEGIN 

START; 

FOR  J  :  =1  TO  NP  DO 
BEGIN 

SETFLAGS  ; 

IF  (IONDIGOl)  OR  (  AMI  J3  0  1.0)  THEN 
BECIN 

JINDEX ; -JINDEX+1 ; 

N  :  -0  ; 

NN : - 0 ; 

1 1  : -0  ; 

IF  JINDEX-INDEX  THEN  VR ITELN( CHR ( 1 2 ) )  ; 

IF  JINDEX-INDEX  THEN  INDEX : - INDEX +5 ; 

LNFD  < 4 )  ; 

WRITELN( FF ,'***•*******•*********••***•*•**•*•***••****•**••'); 
WRITE  (  FF  ,  ' *  AREA  OF  PREDOMINANCE  FOR;  *  >; 

FOR  HH ; - 1  TO  15  DO  WR ITE ( FF , SP C J , HH 3 ) ; 

WRITE  <  FF ,  *  * •  >  ; 

IF  IONDIG-1  THEN  WRITE  (  FF  ,  *  DISSOLVED  SPECIES  DIAGRAM*  >; 
LNFD (  1  )  ; 

WRITELN<  FF ,  •««**********««***«*****•**««*«««***•**«*«**«**«*.). 

FOR  I ; -1  TO  MI  DO 

IF  (LINEC1  , I3-J)  OR  ( L INE C 2 . I  1 - J )  THEN 
IF  < I OND I C  (  )  1  )  OR  ( IONLINC I  1-1 )  THEN 
STABLE2 ; 

STABLE3 ; 

IF  PHMIN< 1 6  THEN 
BEGIN 

STABLE4 ; 

IF  PHMAX  > -  PHMIN  THEN 
BEGIN 

PHI : -PHMIN ; 

PH2 : -PHMAX; 

UPBNDLNS ; 

IF  PASS-0  THEN  FLAG1 : -TRUE ; 

IF  FLAG1-FALSE  THEN 
BEGIN 

IF  PASS-1  THEN  TEMPI; 

WHILE  FLAG1-FALSE  DO  INTPTSHI ; 

END; 

PHI : -PHMIN ; 

PH2 ; -PHMAX ; 

PASS: -0 ; 

LOVBNDLNS ; 

IF  PASS-0  THEN  FLAG3 : -TRUE ; 

IF  FLAG3-FALSE  THEN 
BEGIN 

IF  PASS-1  THEN  TEMP 2 ; 

WHILE  FLAC3-FALSE  DO  INTPTSLO; 

END; 

INTER I MWR I TER ; 

IF  FLAGGER. FALSE  THEN 
BECIN 

STABLES ; 

STABLE* ; 

IF  FLAGS-FALSE  THEN  STABLE7 ; 

IF  (FLAG7-FALSE)  AND  (FLAGS-FALSE)  THEN  STABLES 
IF  ( FLAG7-FALSE )  AND  (FLAGS-FALSE)  THEN  STABLE? 
IF  FLAC7-FALSE  THEN  STABLE  1 0 ; 

IF  FLAC7-FALSE  V  N  STABLE 1 1 ; 

END; 

END; 

END; 

IF  FL ACGER-FALSE  THEN 
BEGIN 

LNFD( 4 ) ; 

WRITE ( FF , '  NO  AREA  OF  PREDOMINANCE'); 

LNFD( 2 ) ; 

END; 

END; 

END; 

LNFD( 4) ; 

END; 


BEGIN 

END. 


<  *  SS+*  ) 


UNIT  INDATA;  INTRINSIC  CODE  18  DATA  19 
INTERFACE 

USES  CONVERT; 

TYPE 

AON*  ARRAY  Cl.. 15,1.. 153  OF  CHAR; 
ASP-ARRAYC 1  .  . 35  ,  1  .  .  1 51  OF  CHAR; 

AO IOkARRAYCI . .10]  OF  REAL; 

A INTI  0* ARRAY  Cl.  .101  OF  INTEGER; 

A3 5* ARRAY C 1 . .353  OF  REAL; 

A I  35  * ARRAY  Cl  .  .351  OF  INTEGER; 

Ml  OX  3  5  * ARRAY  Cl.  .101  OF  A13S; 

VAR 

NAME  : AON ; 

SP  : ASP ; 

T  : REAL ; 

FEHZO  : REAL ; 

NION  : INTEGER; 

F I ON  ; AO 1 0 ; 

A I  ON  : AOl 0 ; 

ZION  : AINTl 0 ; 

B I  ON  ; A I NT 1 0 ; 

OXY  : A  INTI  0 ; 

HYD  ; AINTl 0 ; 

NP  ; INTEGER; 

AM.F  : A35 ; 

Z  : AI 35 ; 

METAL  : AI 35 ; 

OXYCEN; AI3S ; 

HYDRO  : AI 35 ; 

ION  :M10X35; 

FFF  : TEXT ; 

PROCEDURE  GETNONMET; 

PROCEDURE  GETMET; 

PROCEDURE  DATA INPUT; 

PROCEDURE  METAL INPUT; 

PROCEDURE  DATACORRECT; 

PROCEDURE  I ONSOUT 1 ; 

PROCEDURE  IONSOUT2; 

PROCEDURE  NONMETALDI SPLAY ; 
PROCEDURE  DIS1KET; 

PROCEDURE  METALDI SPLAY; 

PROCEDURE  LF(QQ: INTEGER) ; 

PROCEDURE  INCHAR; 

PROCEDURE  INREAL; 

IMPLEMENTATION 
VAR  I; INTEGER; 

Q.M; INTEGER; 

L: INTEGER; 

KK, J: INTEGER; 

K: INTEGER; 

HH: INTEGER; 

XYX : INTEGER; 

PROCEDURE  LF; 

BEGIN 

FOR  L:-l  TO  QQ  DO  VRITELN ; 

END; 

PROCEDURE  DRAWCQQ: INTEGER) ; 

BEGIN 

FOR  L:-l  TO  QQ  DO  WRITE < FFF ,*-•> ; 
WRITELNI FFF ) ; 

END; 


PROCEDURE  INCHAR; 

BEGIN 

FOR  HH : « 1  TO  3  DO  CH3  CHH J  : ■ •  •; 

HH  :  -  1  ; 

WHILE  (NOT  EOLN)  AND  (HH(4>  DO 
BEGIN 

READ ( CH3  C HH I )  ; 

IF  CH3  CHH I »CHR( 8 )  THEN 
BEGIN 

WRITE  <  CHR< 8 ) >  ; 

HH: -HH-2 ; 

END; 

HH: »HH+ 1 
END; 

READLN ; 

CHARTOINT; 

END; 

PROCEDURE  INREAL; 

BEGIN 

FOR  HH : • I  TO  11  DO  CHllCHHl:-1  • 
HH : - 1 ; 

WHILE  (NOT  EOLN)  AND  (HH<12>  DO 
BEGIN 

READ( CHI  1 CHH] )  ; 

IF  CHllCHHl-CHR(S)  THEN 
BEGIN 

WRITE  <  CKR<  6 ) )  ; 

HH: -HH-2 ; 

END; 

HH : *HH+ 1 ; 

END; 

READLN; 

CHARTOREAL ; 

END; 

PROCEDURE  GETNONMET; 

BEGIN 
HH  * ■ 1  * 

WHILE 'NOT  EOLN  DO 
BEGIN 

READ <  NAME  C  K , HH 3 ) ; 

IF  NAME CK . HH] »CHR( 8 )  THEN 
BEGIN 

WR ITE(CHR(8>) ; 

HH: -HH-2 ; 

END; 

HH:  «HH«-1  ; 

END; 

READLN; 

END; 

PROCEDURE  GETMET; 

BEGIN 
HH  :  .  l  ; 

WHILE  NOT  EOLN  DO 
BEGIN 

READ<  SPCM . HH] ) ; 

IF  SP CM >  HH 1 »CHR ( 8 )  THEN 
BEGIN 

WRITE (CHRIS) ) ; 

HH-HH-2 ; 

END; 

HH :  -HH+l  ; 

END; 

READLN; 

END; 


PROCEDURE  DATAINPUT; 

BECIN 

ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 

LF  <  t  >  ; WRITE (  *  ENTER  TEMPERATURE  (IN  DECREES  KELVIN):  •>; 
INREAL : 

END; 

T: -RNUMBER; 

ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 

LF  < 1 )  ; WRITE  <  *  ENTER  CIBBS  FREE  ENERCY  OF  FORMATION'); 

LF ( 1 ) ; WRITE ( '  OF  WATER  AT  ' ,T: 4:2,'  DECREES  (CALS/MOLE):  '); 
INREAL ; 

END  ; 

FEH20 : -RNUMBER; 

ERROR: -TRUE ; 

WHILE  ERROR-TRUE  DO 
BECIN 

LF  < 1 )  ; WRITE (  1  ENTER  THE  NUMBER  OF  NON-METALLIC  SPECIES:  •); 
INCHAR; 

END; 

NION: -NUMBER ; 

IF  NION  <>  0  THEN 

BEGIN 

FOR  I : -1  TO  NION  DO 
BECIN 

K  :  -  I  -4-2  ; 

LF ( 1 )  ; WRITE  < '  NAME  OF  NON-METALLIC  SPECIE  *  ',1:2,':  '); 

CETNONMET; 


END; 

END; 


ERROR -.-TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 

LF ( 1 ) ; WRITE  <  * 
INREAL ; 

END; 

FIONC I J : -RNUMBER; 
ERROR -.-TRUE; 

WHILE  ERROR-TRUE  DO 
BECIN 

LF ( I ) ; WRITE ( ' 
INREAL; 

END; 

AIONC 1 1 : -RNUMBER ; 
ERROR: -TRUE ; 

WHILE  ERROR-TRUE  DO 
BECIN 

LF  < 1 ) .  WRITE ( ' 
INCHAR; 

END; 

ZIONCI1 : -NUMBER; 
ERROR: -TRUE ; 

WHILE  ERROR-TRUE  DO 
BECIN 

LF ( 1 ) ; WRITE  < ' 
INCHAR; 

END; 

BIONC 1 1  .-NUMBER; 
ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BECIN 

LF ( 1 ) ; WRITE ( ' 
INCHAR; 

END; 

OXYCI] : -NUMBER; 

ERROR :  -TRU*. ; 

WHILE  ERROR-TRUE  DO 
BECIN 

LF ( 1 ) ; WRITE ( ' 
INCHAR; 

END; 

HYDCI1 : -NUMBER; 

LF  < 1 > ; 

WRITE LN(CHR< 12 ) ) ; 
END; 


FREE  ENERCY  OF  FORMATION:  '); 


IONIC  ACTIVITY: 


ELECTRICAL  CHARGE:  ’); 


NUMBER  OF  CHARACTERISTIC  ATOMS:  '); 


NUMBER  OF  OXYGEN  ATOMS :  ' > ; 


NUMBER  OF  HYDROGEN  ATOMS: 


PROCEDURE  METAL INPUT; 

BEGIN 

WRITELN<CHR<12> > ; 

ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 

INCHAR?RITE< '  ENTER  NUMBER  OF  METALLIC  SPECIES:  ■>; 

END; 

NP : -NUMBER; 

FOR  M : - 1  TO  NP  DO 
BEGIN 

GET-MET?*1™  '  NAME  °F  METALLIC  SPECIE  •  *  ,  M :  2  ,  •  :  •>; 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 

INREa£^RITE<*  EREE  ENERGY  OF  FORMATION:  •>; 

END; 

FCM] : -RNUMBER ; 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BECIN 

)  ;VRITE<  '  IONIC  ACTIVITY:  •>; 

I NR  CAL  I 
END; 

AMCM1 : -RNUMBER ; 

ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 

INCHAR^RITE<  ELECTRICAL  CHARGE:  ■>; 

END; 

Z  CM]  .-NUMBER; 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BECIN 

INCHAR^RITE< '  NUMBER  OF  METAL  ATOMS:  ■>; 

END; 

METAL  CM! : -NUMBER ; 

ERROR: -TRUE; 

WHILE  ERROR.TRUE  DO 
BEGIN 

r«riiiWRITE< '  NUMBER  OF  OXYGEN  ATOMS:  •>; 

END; 

OXYCENCM] : -NUMBER ; 

ERROR: -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 

INCHARWRITE< '  NUMBER  OF  HYDROCEN  ATOMS:  «>; 

END; 

HYDRO CM] : -NUMBER ; 

IF  NION  >  0  THEN 
BEGIN 

FOR  I : -1  TO  NION  DO 
BEGIN 

ERROR: -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 

ais^'ForiaiSTSKrg  atoms 

INCHAR;  ' 

END; 

10NC I >  Ml : -NUMBER; 

END; 

END 

ELSE  IONC I ,M1 : -0 ; 

LF (  1  )  ; 

WRiTELNf  CHR( t  2  > ) ; 

END; 

END; 


NAME  OF  METALLIC  SPECIE  •  ',M:2, 


FREE  ENERCY  OF  FORMATION:  •>; 


IONIC  ACTIVITY: 


ELECTRICAL  CHARGE: 


NUMBER  OF  METAL  ATOMS:  •>; 


NUMBER  OF  OXYGEN  ATOMS:  •>; 


NUMBER  OF  HYDROGEN  ATOMS: 


PROCEDURE  IONSOUT1; 

BEGIN 

WRITE  <  F  F  F  '  1  1*2  1  !  '); 

FOR  HH  :  ■  l  ' TO  6  *  DO  WRITE  (FFF  ,NAMEC  14-2  ,  HH)  )  ; 

WRITE  (FFF,  ’  I  '  ,  F  I  ON  111:10:1,  ' 

WRITELNI  FFF  ,  A I  ONC  11:2:4,'  !  '  ,  Z  I ONC  I  I  :  4  ,  '  I*); 

END  ; 

PROCEDURE  IONSOUT2 ; 

BEGIN 

WR ITE  <  FFF ,  *  *,1:2.  •  !  *  , B I ONC 1 1 : 3 ,  * 

WRITELN<  FFF  ,  1  ',OXYCU:2,'  !  ',HYDCIJ:2,' 

END; 

PROCEDURE  NONMETALDI SPLAY; 

BEGIN 

WRITELNI CHR< 1 2 ) ) ; 

WRITELN( FFF ) ; 

WRITELN< FFF , ***************************** ;3?>; 
WRITELN( FFF ,  *  *  TEMPERATURE  -  ':29,T:S:2,'  *'); 

WRITELN( FFF , •******«******************«*': 39); 
WRITELN(FFF) ; 

IF  NION-O  THEN 
BEGIN 

WRITELN(FFF) ; WR ITELN ( FF F ) ; 

WRITE  <  F  FF ,  ' *  *  *  NO  NON-METALLIC  IONIC  SPECIES 
WRITE!FFF ,  *  ARE  CONSIDERED  ***•>; 

N I ON : « 1 ; 

END 

ELSE 

BEGIN 

WRITELN(FFF) ; 

DRAW (S3) ; 

WRITE  <  FFF ,  '  ION  !  IONIC  !  FREE  ENERGY  •); 
WRITELN!FFF, ’ !  IONIC  !  I’>; 

WRITE! FFF, 'NUMBER  I  SPECIES  1  OF  FORMATION  '); 
WRITELN<FFF. ' I  ACTIVITY  I  CHARCE  I'); 

DRAW! S3) ; 

FOR  l:«l  TO  NION  DO 
I ONSOUT1 ; 

FOR  I: -l  TO  3  DO  WRITELN!FFF> ; 

DRAW ! 4  5 )  ; 

WRITELN!FFF,  '  IO.I  J  CHARACTERISTIC  1  OXYGEN 
WRITELN! FFF, 'NUMBER  !  ATOMS  I  ATOMS 

DRAW! 45); 

FOR  I ; -1  TO  NION  DO 
IONSOUT2 ; 

END; 

END; 

PROCEDURE  DI St MET ; 

BEGIN 

WRITE! FFF,  '  ' )  ; 

FOR  HH : ■ 1  TO  IS  DO  WRITE  I FFF , SPC I , HKJ > ; 

WRITE!FFF ,  '  I  '  , METAL C I  1 : 2 .  '  I  '  , OXYGENC 1 1 : 2 > ; 
WRITE! FFF , '  \  ' , HYDRO C l 3 : 2 , '  I'); 

WRITELN! FFF) ; 

END; 

PROCEDURE  DIS2MET; 

BEGIN 

WRITE ! FFF , '  '); 

FOR  HH ; ■ I  TO  IS  DO  WRITE  I FFF , SPt I ,HHJ ) ; 

WRITE! FFF , ' 1 ' > ; 

IF  Q-I  THEN 

FOR  J:-l  TO  S  DO 
IF  J <  «N I ON  THEN 

WRITE  !  FFF  ,  '  '  ,  I  ON  C «? ,  1 1  ;  2  ,  '  I'); 

IF  Q>2  THEN 

FOR  J TO  10  DO 
IF  J ! aN I ON  THEN 

WRITE ! FFF .  '  ' , I ON  C  J , 1 1 : 2 , '  I'); 

IF  Q-3  THEN 

FOR  J:-!l  TO  IS  DO 
IF  J ! aNION  THEN 

WRITE ! FFF , *  ' , IONC J. 1 1 : 2 , '  I'); 

IF  Q-4  THEN 

FOR  J:-l*  TO  20  DO 
IF  J  !  *>W  I  ON  THEN 

WRITE ! FFF , '  ' , IONC J , I J : 2 , '  I'); 

WRITELN! FFF) ; 

END; 


!  '  >  ; 

!  '  >  ; 


)  ; 


HYDROGEN  l ' ) 
ATOMS  ! ' ) 


PROCEDURE  DIS3MET; 

BEGIN 

WRITE  (FFF,  1  •  )  ; 

IF  Q>1  THEN 


WR ITE ( FFF , ' ! 

I  ONI  1 

1 0N2  < 

I  ON  3  ! 

I  ON  4  ! 

IONS  1 

) 

IF  Q-2  THEN 
WRITE ( FFF , ' ! 

I0N6  ! 

I  ON  7  1 

IONS  ! 

I  ON?  ! 

IONI 0  t 

) 

IF  0-3  THEN 
WRITE ( FFF , ' ! 

I  ON  1  1  ! 

I ON l 2 i 

I  ON 1 3 ! 

I ONI  4 1 

IONI S ! 

) 

IF  0-4  THEN 
WRITE (FFF , ' ! 

I ONI  6 < 

I  ONI  7 1 

I ON 1 8 ! 

IONI 9 1 

ION20 ! 

) 

WRITELN! FFF  >  ; 

DRAW  (S3)  ; 

END; 

PROCEDURE  METALDISPLAY ; 

BEGIN 

WRITELN! CHR< 12  >  )  ; 

WRITELN(FFF) ; WR ITELN < F F F > ; 

DRAW (33) ; 

WRITE  LN  ( FFF,'  !  FREE  ENERGY  !  SPECIES  I 

WRITELN! FFF , 'METALLIC  SPECIES  !  OF  FORMATION  I  ACTIVITY  !  CHARGE 
DRAW! S3) ; 

FOR  I:-l  TO  NP  DO 
BECIN 

WRITE ( FFF  '  *) • 

FOR  HH :  - 1  ' TO  1 5 ' DO  WR ITE ( FFF . SP C I , HH ) > ; 

WRITE (FFF,  ‘  I  ‘  ,FC 11 : 11  : l ,  ‘  !  ’  . AMt I  3  :  2  :  4 >  ; 

WRITELNtFFF, *  I  ',ZCI3:2,'  !‘>; 

END; 

FOR  I:-l  TO  3  DO  WRITELN(FFF) ; 

DRAW! 44) ; 

WRITELN! FFF, ‘  !  METAL  I  OXYGEN  1  HYDROGEN  !‘>; 

WRITELN(FFF . 'METALLIC  SPECIES  I  ATOMS  I  ATOMS  I  ATOMS  !*>; 
DRAW<  4  6 )  ; 

FOR  I:-l  TO  NP  DO  DISIMET; 

FOR  I:-l  TO  3  DO  WRITELN! FFF > ; 

L ; -TRUNC! (NION-l ) /5 . 0  >♦!  ; 

FOR  Q; -1  TO  L  DO 
BEGIN 

WRITELN(FFF) ; 

DRAW (53); 

WRITE  (FFF  1  •  )•• 

WRITELN! FFF , *  1  NUMBER  OF  CHARACTERISTIC  ATOMS  !‘>; 

WRITEIFFF, 'METALLIC  SPECIES  I  OF  EACH  IONIC  SPECIES'); 

WRITELN! FFF, '  I ' ) ; 

DIS3MET; 

FOR  I:-l  TO  NP  DO  DIS2MET; 

FOR  I  :  •  l  TO  3  DO  WR ITELN! FFF > ; 

END; 

END; 

PROCEDURE  NMETC7 ; 

BEGIN 

WRITE!'  MODIFICATION  OF  NUMBER  OF  HYDROGEN  ATOMS  FOR  '); 

FOR  K:-l  TO  IS  DO  WRITE (NAME C KK , K1 ) ; 

LF ( 3  )  ; 

WRITELNt  '  CURRENT  VALUE  IS  ',KYDCI3:3); 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  (  1  )  ; 

WRITE!'  DESIRED  VALUE: 

INCHAR; 

END; 

KYDC 1 1 : -NUMBER; 

END; 


PROCEDURE  NMETC6 ; 

BEGIN 

WR ITE ( ‘  MODIFICATION  OF  NUMBER  OF  OXYGEN  ATOMS  FOR  •>; 

FOR  K:«l  TO  15  DO  WR I TE < NAME C KK , K 1 > ; 

LF  <  3  )  ; 

WRITELN!  *  CURRENT  VALUE  IS  ',OXYCI3:3>; 

ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  <  1  )  ; 

WR ITE  <  '  DESIRED  VALUE:  '); 

INCHAR; 

END; 

OX  Y  C I  3  : -NUMBER ; 

END; 

PROCEDURE  NMETC5 ; 

BEGIN 

WRITE  < '  MODIFICATION'  OF  NUMBER  OF  CHARACTERISTIC  ATOMS') 

FOR  K : - 1  TO  15  DO  WR ITE ( NAME C KK , KI > ; 

LF  <  3  )  ; 

WRITELN<  '  CURRENT  VALUE  IS  '  , BIONC I  3 : 3 ) ; 

ERROR: -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  < 1  )  ; 

WRITE! '  DESIRED  VALUE:  '); 

INCHAR; 

END; 

BIONC I  3 : -NUMBER ; 

END; 

PROCEDURE  NMETC4 ; 

BEGIN 

WRITE! '  MODIFICATION  OF  CHARGE  FOR  *>; 

FOR  K:  -1  TO  15  DO  WRITE ! NAME C KK , K 3 > ; 

LF  !  3  )  ; 

WRITELN! '  CURRENT  VALUE  IS  ' , ZIONC 11:4); 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF!  I  >  ; 

WRITE!'  DESIRED  VALUE:  '); 

INCHAR; 

END; 

ZIONC I  3 : -NUMBER ; 

END; 

PROCEDURE  NMETC3 ; 

BEGIN 

WRITE!'  MODIFICATION  OF  IONIC  ACTIVITY  FOR  '); 

FOR  K : - l  TO  IS  DO  WRITE !NAME C KK , K3 ) ; 

LF  !  3  )  ; 

WRITELN!  '  CURRENT  VALUE  IS  *  , AIONC 13:2:6); 

ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF!  1)  ; 

WRITE!'  DESIRED  VALUE:  '); 

INREAL; 

END  ; 

AIONC I  3 : -RNUMBER ; 

END; 

PROCEDURE  NMETCZ ; 

BEGIN 

WRITE!'  MODIFICATION  OF  FREE  ENERGY  FOR  '); 

FOR  K : ■ 1  TO  IS  DO  WR ITE !NAME C KK , Kl ) ; 

LF!  3)  ; 

WRITELN! '  CURRENT  VALUE  IS  ' , FIONC 13:10:2); 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF!  1  > ; 

WRITE!'  DESIRED  VALUE:  '); 

INREAL; 

END; 

FIONC I  3 : -RNUMBER; 

END; 


PROCEDURE  NMETCl ; 

BEGIN 

WR ITELN <  1  MODIFICATION  OF  SPECIES  NAME  *  > ; 

LF  !  2  )  ; 

WRITE! *  CURRENT  NAME  IS  *>; 

FOR  K : » 1  TO  15  DO  WR I TE ! NAME C KK , K 1 > ; 

LF ( 2  >  ; 

WRITE! '  DESIRED  NAME:  •>; 

FOR  K: >1  TO  15  DO  NAME  II . K]  :  -  "  •; 

K. :  «KK  ; 

GETNONMET; 

END; 

PROCEDURE  NONMETCOR; 

VAR  FLAG2 : BOOLEAN; 

BEGIN 

WR ITELN! CHR ! 1 2  > )  ; 

WR ITELN!  '  MODIFICATION  OF  NON-METALLIC  SPECIES  DATA 1  > 

ERROR: -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  !  3  )  ; 

WRITE!*  ENTER  IONIC  SPECIE  NUMBER:  *>; 

INCHAR; 

IF  (NUMBER! 1 )  OR  ( NUMB  E  R  >  N I  ON )  THEN  ERROR : -TRUE ; 

END; 

I : -NUMBER; 

KK  :  -  I  4-2  ; 

REPEAT 

FLAG2 : -FALSE ; 

LF ( 2  )  ; 

WRITELN! '  PARAMETER  TO  BE  MODIFIED  -  *); 

LF ! 1 >  ; 

WRITELN! *  1  -  ION  NAME'); 

WRITELN! *  2  -  FREE  ENERGY  OF  FORMATION* >; 

WRITELN! '  3  -  ACTIVITY* >; 

WRITELN!*  4  -  CHARGE*); 

WRITELN!'  5  -  NUMBER  OF  CHARACTERISTIC  ATOMS*); 

WRITELN! *  6  -  NUMBER  OF  OXYGEN  ATOMS*); 

WRITELN! *  7  -  NUMBER  OF  HYDROGEN  ATOMS*); 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF! 1 ) ; 

WRITE!'  ENTER  OPTION:  *); 

INCKAR; 

END; 

J : -NUMBER; 

IF  (J!l>  OR  ( J>  7 )  THEN  FL AG2 : -TRUE ; 

UNTIL  FLAG2-FALSE ; 

WRITELN! CHR! 12) ) ; 

CASE  J  OF 

1:  NMETCl; 

2 :  NMETC2 ; 

3:  NMETC3 ; 

4 :  NMETC4 ; 

5 :  NMETC5 ; 
i :  NMETC4 ; 

7 :  NMETC7 ; 

ENO; 

END; 


PROCEDURE  METC8 ; 

BEC  IN 

WRITE! '  MODIFICATION  OF  *  >; 

FOR  K:-l  TO  15  DO  WR ITE ! SP C I , K 3 > ; 

ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  <  3  >  ; 

WRITE! 1  CHANGE  VALUE  FOR  WHAT  ION  NUMBER?:  •>; 

INCHAR; 

IF  ! NUMBER  < 1 )  OR  < NUMB ER > N I  ON )  THEN  ERROR : -TRUE ; 

END; 

K: -NUMBER; 

LF! 2  >  ; 

WR ITELN!  ‘  CURRENT  VALUE  IS  *  ,  I ONC K , I  3  : 2 >  ; 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  !  1  >  ; 

WRITE! '  DESIRED  VALUE:  '); 

INCHAR; 

END; 

IONCK,  13  -.-NUMBER; 

END; 

PROCEDURE  METC7 ; 

BEGIN 

WRITE!'  MODIFICATION  OF  NUMBER  OF  HYDROGEN  ATOMS  FOR 

FOR  K  :  - 1  TO  15  DO  WR ITE ! SP Cl , K 3 > ; 

LF  !  3  )  ; 

WRITELN!  '  CURRENT  NAME  IS  '  . HYDRO C I  3  : 3 >  ; 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF!  1  )  ; 

WRITE!'  DESIRED  VALUE:  '); 

INCHAR: 

END; 

HYDROC I  3  : -NUMBER ; 

END; 

PROCEDURE  METC8 ; 

BEGIN 

WRITE!'  MODIFICATION  OF  NUMBER  OF  OXYGEN  ATOMS  FOR  ') 

FOR  K : - 1  TO  15  DO  WRITE  I SPC I , K3 ) ; 

LF  !  3  >  ; 

WRITELN!'  CURRENT  VALUE  IS  ' , OX YGENC 13:3); 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF!  1  )  ; 

WRITE!'  DESIRED  VALUE:  '); 

INCHAR; 

END; 

OXYGENC I  3 : -NUMBER ; 

END; 

PROCEDURE  METC5 ; 

BEGIN 

WRITE!'  MODIFICATION  OF  NUMBER  OF  METAL  ATOMS  FOR  '>; 

FOR  K : -1  TO  15  DO  WRITE ! SPC I , K3 )  ; 

LF  !  3  )  ; 

WRITELN!'  CURRENT  VALUE  IS  '  , METAL C I  3  . 3 > ; 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  !  1  >  ; 

WRITE!'  DESIRED  VALUE:  •>; 

INCHAR; 

END; 

METAL I  I J : -NUMBER; 

END; 


PROCEDURE  METC4 ; 

BEGIN 

WR ITE  <  '  MODIFICATION  OF  CHARCE  FOR  '); 

FOR  K  :  ■  I  TO  15  DO  WRITE ! SPC I . KI > ; 

LF  (  3  >  ; 

WRITELN! '  CURRENT  VALUE  IS  ',ZCI3:2>; 

ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  !  1  >  ; 

WRITE  < '  DESIRED  VALUE:  •>; 

INCHAR; 

END; 

ZC  I  3  : -NUMBER; 

END; 

PROCEDURE  METC3 ; 

BECIN 

WRITE! •  MODIFICATION  OF  ACTIVITY  FOR  '  > 

FOR  K:-l  TO  15  DO  WR ITE ! SP C I , K 3 >  ; 

LF  <  3  >  ; 

WR ITELN< '  CURRENT  VALUE  IS  *  , AMI  I  3 : 2 : 6 )  ; 
ERROR; -TRUE; 

WHILE  ERROR-TRUE  DO 
BECIN 
LF!  1  >  ; 

WRITE! '  DESIRED  VALUE: 

INREAL ; 

END; 

AMC I  1 : -RNUMBER ; 

END; 

PROCEDURE  METC2 ; 

BECIN 

WRITE! •  MODIFICATION  OF  FREE  ENERCY  FOR 

FOR  K : - 1  TO  15  DO  WRITE < SPC I , Kl > ; 

LF  !  3  >  ; 

WRITELN!  ‘  CURRENT  VALUE  IS  '  .  F  C  I  3  :  1  0  :  2  >  ; 
ERROR; -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  < 1  )  ; 

WRITE!’  DESIRED  VALUE:  ’); 

INREAL; 

END; 

FC 1 1 : -RNUMBER; 

END; 

PROCEDURE  METC1 ; 

BEGIN 

WRITELN! '  MODIFICATION  OF  SPECIE  NAME') 

LF! 3  >  ; 

WRITE!'  CURRENT  NAME  IS  '); 

FOR  K :  - 1  TO  15  DO  WRITE ! SPC I , K I > ; 

LF!  2)  ; 

WRITE!'  DESIRED  VALUE:  •>; 

FOR  K : - 1  TO  IS  DO  SPCI.JC3:-' 

M:  -I  ; 

GETMET; 

END; 


•  « 


PROCEDURE  METCOR; 

VAR  FLAC 1 : BOOLEAN; 

BECIN 

ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BECIN 
LF  <  1  >  ; 

WRITE ( '  ENTER  METALLIC  SPECIE  NUMBER:  •>; 

I NCH AR ‘ 

IF  ( NUMBER < 1 '  OR  <NUMBER>ND  THEN  ERROR : -TRUE ; 

END; 

I . -NUMBER ; 

REPEAT 

FLAC 1 : -FALSE ; 

LF  (  2  )  ; 

WR ITELN <  •  PARAMETER  TO  BE  MODIFIED  -•>; 

LF  <  1  >  ; 

WR ITELN< ‘  1  -  SPECIES  NAME’); 

WRITELNI *  2  -  FREE  ENERGY  OF  FORMATION1 >; 

WRITELN< 1  3  -  ACTIVITY1); 

WRITELNI 1  4  -  CHAFGE 1 ) ; 

WRITELNI 1  5  -  NUMBER  OF  METAL  ATOMS1); 

WRITELNI1  6  -  NUMBER  OF  OXYCEN  ATOMS1); 

WRITELNI 1  7  -  NUMBER  OF  HYDROGEN  ATOMS1); 

WRITELNI*  8  -  NUMBER  OF  CHARACTERISTIC  ATOMS1); 

ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BECIN 
LF  I  1 )  ; 

WRITE!1  ENTER  OPTION:  ‘>; 

INCHAR; 

END; 

J : -NUMBER ; 

IF  IJ<1)  OR  I  J ) 8 )  THEN  FLAC 1  : -TRUE ; 

UNTIL  FLAC1 -FALSE ; 

WRITELNI CHR! 12) ) . 

CASE  J  OF 

1:  METC1; 

2:  METC2 ; 

3:  METC3; 

4 :  METC4 ; 

5 :  METC5 ; 

6 :  METC4 ; 

7:  METC7 ; 

8 :  METC8 ; 

END; 

END; 

PROCEDURE  TEMPCOR ; 

BECIN 

WRITELNI*  CURRENT  TEMPERATURE  VALUE  IS  *,T:4:2,‘  DECREES  KELVIN1) 
ERROR: -TRUE; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  I  I  )  ; 

WRITE!1  DESIRED  VALUE:  *>; 

INREAL; 

END; 

T: -RNUMBER; 

END; 


PROCEDURE  FECOR ; 

BEGIN 

WRITELNt '  CURRENT  FREE  ENERGY  IS  ' , FEH20 : 8 : 2 , '  CALS/MOLE 
ERROR : -TRUE ; 

WHILE  ERROR-TRUE  DO 
BEGIN 
LF  <  1  )  ; 

WRITE! ‘  DESIRED  VALUE:  •>; 

INREAL ; 

END  i 

FEH20: -RNUMBER; 

END; 

PROCEDURE  DATACORRECT; 

VAR  OPT: CHAR; 

BEGIN 

WRITELNt CHR< 12) )  ; 

REPEAT 

WR ITELN ( CHR  < 12)); 

LFt  1  >  ; 

WRITELNt '  MODIFY  WHICH  VALUE?1); 

LF(  1  >  ; 

WRITELNt1  1  -  TEMPERATURE'); 

WRITELNt 1  2  -  FREE  ENERGY  OF  H201); 

WRITELNt 1  3  -  NON-METALLIC  ION  DATA1); 

WRITELNt1  4  -  METALLIC  SPECIES  DATA1); 

WRITELNt1  5  -  ALL  DATA  IS  NOW  CORRECT1); 

LFt  1  )  ; 

WRITEt1  ENTER  OPTION:  ’); 

RE ADLNt  OPT  >  ; 

WR ITELNt  CHRt 12 )  )  ; 

CASE  OPT  OF 

1 1 1 :  TEMPCOR; 

1 2 1 :  FECOR ; 

'3':  NONMETCOR; 

1 4 1 :  METCOR; 

•S':  OPT : -  1 N 1  ; 

END; 

LF  t  2  )  ; 

UNTIL  OPT- ‘N1 ; 

END; 


BEGIN 

END. 


USER'S  GUIDE  FOR  THE 
INTERACTIVE  PH  -  POTENTIAL  PROGRAM 


This  introductory  guide  should  be  read  prior  to  operation 


of  the  program 


When  all  else  fails.  push  the  reset  button 

When  the  reset  button  fails*  pull  the  plug 


Eyinhard  Muth 


INTRODUCTION 


The  microcomputer  version  of  the  PH  -  POTENTIAL  PROGRAM  is 
easy  to  use.  even  for  those  with  little  or  no  computer  experience. 
This  guide  explains,  in  simple  terms,  how  to  get  started  and  what 
the  operator  must  do  to  use  the  program  effectively. 


GETTING  STARTED 

Getting  started  is  easy.  Simply  place  the  diskettes  in  the  proper 
disk  drives  and  turn  on  the  power.  The  program  will  begin  execu¬ 
tion  automatical ly .  The  following  explains  the  process  in  detail. 

1  -  Open  disk  drive  #1.  Insert  the  diskette  labeled  POURBAIX  1 

into  disk  drive  #1.  (For  the  beginner,  hold  the  diskette 

in  either  hand,  palm  up.  and  with  your  thumb  on  the  label, 

gently  slide  in  the  diskette  and  close  the  disk  drive  door.) 

2  -  Repeat  step  1  for  disk  drive  #2  and  the  diskette  labeled 

POURBAIX  2. 

NOTE:  For  single  drive  systems,  the  program  can  be  ex¬ 
ecuted  with  only  POURBAIX  1. 

3  -  Turn  on  the  power  for  the  computer.  If  the  power  was 

already  on  turn  the  computer  off  and  then  on  again.  The 
disk  drives  will  make  whirring  noises  as  the  computer  pre¬ 
pares  for  operation.  In  just  a  few  seconds,  the  program 
will  begin  execution  automatically  and  prompt  for  inputs. 


ENTERING  THE  DATA 


The  program  sequentially  prompts  for  the  required  data-  Simply 
type  in  the  name  or  number  as  asked.  If  the  computer  responds 
with  "IMPROPER  NUMBER"  see  Appendix  1  for  possible  causes. 

Data  to  be  Entered  (R  =  Real.  I  =  Integer.  C  =  Character. ) 

1  -  TEMPERATURE  (degrees  Kelvin).  <R> 

2  -  GIBB 'S  FREE  ENERGY  OF  FORMATION  OF  H20  AT  THE  TEMPER¬ 

ATURE  BEING  CONSIDERED  (in  calories  per  mole).  (R> 

3  -  NUMBER  OF  NON-MET ALL I C  IONS  CONSIDERED.  ( I > 

For  each  non  -metallic  ionic  species  considered: 

4A  -  NAME  OF  NQN-METALLIC  ION.  (C> 

4B  -  FREE  ENERGY  OF  FORMATION  (in  calories  per  mole).  (R> 
4C  -  THERMODYNAMIC  ACTIVITY.  (R> 

4D  -  ELECTRICAL  CHARGE.  (  I  > 

4E  -  NUMBER  OF  CHARACTERISTIC  ATOMS  (atoms  other  than 
oxygen  or  hydrogen).  (!) 

4F  -  NUMBER  OF  OXYGEN  ATOMS  IN  THE  IONIC  SPECIES.  ( I > 

4C  -  NUMBER  OF  HYDROGEN  ATOMS  IN  THE  IONIC  SPECIES.  < I > 


3  -  NUMBER  OF  METALLIC  SPECIES  CONSIDERED.  (I> 


For  each  metallic  -  .  cies  considered: 

6A  -  NAMH  OF  ME  I  ALL  I C  SPECIES.  <C) 

6B  -  GIBB'S  FREE  ENERGY  OF  FORMATION  (calories  per  mole).  <R> 
6C  -  THERMODYNAMIC  ACTIVITY.  <R> 

6D  -  ELECTRICAL  CHARGE.  < I > 

6E  -  NUMBER  OF  METAL  ATOMS  IN  THE  SPECIES.  <I> 

6F  -  NUMBER  OF  OXYGEN  ATOMS.  < I > 

6G  -  NUMBER  OF  HYDROGEN  ATOMS.  < I > 

6H  -  NUMBER  OF  CHARACTERISTIC  ATOMS  OF  IONIC  SPECIES  #1, 

2,  3.  etc.  (I> 

PROGRAM  EXECUTION 

Execution  of  the  program  begins  with  the  last  piece  of  data 
entered.  The  following  message  will  appear  on  the  screen: 

MODIFY  THE  DATA  ?  CY/N3: 

Type  N  if  the  data  was  entered  correctly  or  if  you  want  to 
see  the  data  before  modification. 

Typing  Y  will  start  the  data  modification  routine.  The  computer 
will  then  display  the  following  menu: 

MODIFY  WHICH  VALUE  ? 

1  -  TEMPERATURE 

2  -  FREE  ENF:RCY  OF  H20 

3  -  MON-METALLIC  ION  DATA 

4  -  METALLIC  SPECIES  DATA 

3  -  ALL  DATA  IS  NOW  CORRECT 


ENTER  OPT  ION 


For  options  1  and  2.  the  old  values  are  displayed  along 
with  a  prompt  for  the  desired  value.  For  options  3  and  4. 
the  computer  prompts  for  the  parameter  to  be  modified  before 
displaying  the  old  value  and  awaiting  input  of  the  new  value. 
Once  a  value  has  been  modified,  the  "MODIFY  WHICH  VALUE  ?" 
menu  is  again  displayed.  Option  S  will  end  the  data  modifi¬ 
cation  routine. 

The  computer  will  next  prompt  with: 

DISPLAY  TADULATED  DATA  ?  CY/NJ: 

Type  N  if  you  do  not  want  to  see  the  data  displayed  in  tabular 
form. 

Typing  Y  will  start  the  data  display  routine.  The  computer 
prompts  with: 

DATA  DISPLAY  OPTIONS 

P  -  PRINTER 
T  -  TERMINAL 

ENTER  OPTION: 

The  data  is  then  displayed  as  specified. 

APPLE  NOTE:  To  temporarily  halt  execution  of  the  program  type 
<  CTRL  >  S  i  that  is  type  S  while  depressing  the  CTRL  key. 

This  freezes  the  screen  display  until  (CTRL)  S  is  typed  a¬ 


gain 


The  computer  next  prompts  with: 


EQUATION  LISTING  OPTIONS 

P  -  PRINTER 
T  -  TERMINAL 
ENTER  OPTION: 

The  equations  will  be  listed  as  specified. 

The  computer  displays  the  message  "EXECUTION  COMPLETE...."  to 
signal  the  end  of  the  run  and  prompts  with: 

EXECUTE  AGAIN  ?  EY/N3: 

Type  N  to  quit. 

Type  Y  and  the  computer  asks: 

REPEAT  EXECUTION  WITH  SAME  DATA  ?  CY/N3: 

Type  N  and  the  computer  asks: 

MODIFY  CURRENT  DATA  SET  ?  CY/NJ : 

Type  Y  to  initiate  the  data  modification  routine. 

Type  N  and  the  system  reinitializes  and  prompts  for  new 


inputs 


APPENDIX  1 

REASONS  FOR  "IMPROPER  DATA"  MESSAGE 


1  -  Entering  an  alphabetic  character  for  a  numeric  variable. 

Retype  the  number  correctly. 

2  -  A  decimal  point  (. )  for  an  integer  variable. 

Retype  the  number  omitting  the  decimal  point. 

3  -  A  plus  sign  <  +  )  for  electric  charge. 

Retype  the  number  omitting  the  plus  sign. 

Note:  Negative  values  do  require  the  minus  sign. 


PROGRAMMING  NOTES 


The  interactive  version  of  the  PH  -  POTENTIAL  PROGRAM  is  largely 
a  translation  of  the  WATFIV  version  that  makes  use  of  UCSD  - 
APPLE  PASCAL  features.  Despite  the  limitations  imposed  on  the  size 
of  procedures  and  the  limitations  imposed  by  the  stack,  the  program 
is  designed  in  modules  to  facilitate  modification  and/or  extensions. 
The  program  is  intended  to  be  imp  1 ementab le  on  any  microcomputer 
that  supports  the  UCSD  operating  system. 

The  current  version  makes  extensive  use  of  Intrinsic  Units,  sep¬ 
arately  compiled  program  segments  that  have  been  installed  into 
the  SYSTEM.  LIBRARY.  These  Units  are  suiapped  into  active  memory 
only  uihen  a  call  is  made  to  a  procedure  found  in  the  Unit.  A 
brief  description  of  each  of  the  Units  follows. 

GET  VAL 

This  Unit  converts  an  integer.  CCj3»  into  a  corresponding  character 
array.  SAVE.  For  example,  a  CCj3  value  of  123  is  converted  to  the 
array  of  values  SAVEC 1 3= ' 1 '.  SAVEC23='2'«  SAVEC33='3'. 

CONVERT 

This  Unit  converts  character  arrays  to  integer  or  real  values 
through  the  procedures  CHARTOINT  and  CHARTOREAL. 


INDATA 

This  Unit  contains  the  data  entry 
All  inputs  are  read  as  characters 
type  mismatches.  Integer  and  real 
VERT  Unit. 


and  data  modification  routines, 
to  prevent  crashing  by  variable 
values  are  extracted  by  the  CON¬ 


STABLE 

This  Unit  calculates  the  areas  of  predominance  for  the  species 
being  considered. 


METALLIC  SPECIES 


NUMBER  OF  CHARACTERISTIC  ATOMS  ! 

OF  EACH  IONIC  SPECIES  1 

t  ONI  !  I ON2  I  I ON3  I  ION4  I  IONS  ! 


FE 

FE304 
F  E  203 
HFE02- 
FE0H+4 
FE  <  OH  >  2 ♦ 
FECL2  < A0> 
F  ECL  ♦  ♦ 
FES<  SOL  > 

F  ES2 ( SOL ) 


NUMBER 


EQUILIBRIUM 


BALANCED  CHEMICAL  EQUATION 


i  i 


t 

1 

E  = 

-0  0864 

- 

0  0  5  9  2  PH 

I 

1 

2 

1 

1 

1 

E  = 

-0.0530 

- 

0  0592PH 

I 

9 

i 

3 

1 

1 

1 

E  = 

0  3229 

- 

0  0  8  8  7PH 

i 

i 

4 

1 

1 

1 

E  = 

-0  09  1  4 

- 

0  0  1  9  7  PH 

i 

i 

i 

5 

1 

1 

1 

E  ** 

0  0  0  7  1 

- 

0  .  0  3  9  4  PH 

i 

i 

i 

6 

1 

1 

1 

E  = 

-0.5270 

i 

t 

! 

7 

1 

1 

1 

E  = 

-01424 

i 

i 

9 

a 

1 

1 

1 

E  - 

-05532 

- 

0  .  0  2  9  6  PH 

9 

; 

i 

9 

1 

1 

1 

E  = 

-0.4654 

- 

0  .  029  6PH 

? 

! 

i 

10 

1 

I 

! 

E  = 

0.0878 

- 

0 . 0532PH 

i 

9 

i 

1  1 

» 

i 

i 

E  = 

0.2144 

- 

0 . 0592PH 

i 

9 

9 

1  2 

! 

1 

I 

E  “ 

-13144 

♦ 

0 . 0296PH 

1 

9 

9 

1  3 

i 

l 

l 

E  = 

-0.1318 

♦ 

0  .  2958PH 

i 

9 

9 

14 

t 

I 

l 

E  ■= 

0.7355 

♦ 

0  .  1  1  8 3  PH 

i 

i 

i 

13 

l 

l 

1 

E  = 

1  .2333 

- 

0  .  2367PH 

9 

9 

{ 

16 

1 

1 

1 

E  * 

-0.3902 

♦ 

0  .  4  7  3  3  PH 

t 

! 

I 

17 

l 

l 

1 

E  = 

1  .3141 

- 

0 . 1 479PH 

i 

i 

i 

18 

l 

i 

i 

£» 

-12235 

♦ 

0 . 0296PH 

9 

i 

9 

17 

« 

i 

i 

E  ■ 

0.1511 

- 

0  03  1  1  PH 

i 

! 

! 

20 

i 

i 

! 

E  ■ 

-0.8048 

9 

9 

9 

2  1 

i 

1 

i 

PH- 

0 

.  98 

i 

i 

! 

22 

i 

i 

PH- 

-3 

.  05 

i 

! 

I 


3  *  F  E  ♦  4  *  H  2  0  * 

FE304  ♦  8  *H*  ♦  8*E- 

2  *  F  E  ♦  3  *  H2  O  - 

FE203  4  4  * H ♦  ♦  6  *  E  - 


FE 

♦  2  *  H20  - 

HFE02-  ♦ 

3*H+  +  2  *  E  - 

FE 

♦  H20  - 

F  EQH«-  4  ♦ 

H+  +  3  *  E  - 

FE 

♦  2  *  H2  0 

F  E  (  OH  >24 

+  2*H4  4  3  *  E - 

FE 

♦  2  *  CL- 
FECL2 ( AQ) 

4  2  *  E  - 

FE 

4-  CL- 
FECL+4  ♦ 

3  *  E  - 

FE 

4-  HS- 
FES(SOL) 

4  H  4  4  2  *  E  - 

FE 

4-  2  *  HS- 

FES2 (SOL) 

4  2*H4  4  4*E 

FE 

4  4*H20  ♦ 

4-  9*H4- 

HS- 

4  10*E- 

2  *  FE 304  4  H20 

3  *  F  E  2  03  4  2  *H+  ♦  2*E- 

3*HFE02-  -  2  *  H20  - 

FE304  -  H+  ♦  2  *  E  - 

FE304  -  H20  - 

3  *  FE0H4  4  -  5  *H*  4  E- 

FE304  ♦  2  * H2 0  . 

3  *  FE  <  OH  >24  -  2  *H4  4  E- 

3*FECL2< AQ  >  ♦  4*H20  -  6*CL- 

FE304  4  8*H4  4  2  *  E- 

F E 304  -  4  *H20  4  3*CL- 

3  * FECL 4 ♦  -  8«H4  ♦  E- 

3  *  F ES  <  SQL )  ♦  4  *H20  -  3*HS- 

FE304  4  5*H4  ♦  2  *  E- 

F  E  304  -  4  *  H2  O  ♦  6*HS-  - 

3  *  FES2 ( SOL )  -  2*K+  ♦  4«E 

F E 304  4  8  *H20  ♦  3«HS-  - 

3*  ♦  19*H4  ♦  2  2  *  E- 

2  *HF  E02 -  -  K20  - 

FE203  ♦  2  <  E - 

2  «  FEOH  +  ♦  4  H20  . 

FE203  4  4  *  H ♦ 

2  *  F E ( OH >24  -  H20  - 

FE203  ♦  2«H4 


2  3 

1 

1 

E  * 

0  8  7  5  0 

-  0  1  7  7  5  PH 

1 

1 

2*FECL2(6Q)  ♦  3*H20  -  4*CL- 

EE203  ♦  6  *  H ♦  ♦  2  *  E - 

2  4 

I 

1 

1 

PH« 

1.51 

1 

l 

| 

2*FECL++  ♦  3»H20  -  2»CL- 

FE203  ♦  6  *  H ♦ 

25 

1 

1 

E  * 

0  7  476 

-  0  1  1  8  3  PH 

1 

1 

2  *  F  ES ( SOL )  ♦  3  *  H  20  -  2*HS- 

FE203  ♦  4  *  H ♦  4  2«E- 

26 

1 

1 

1 

E  = 

-  1  7028 

4  0  0  5  7  2  PH 

f 

1 

FE203  -  3  *H2  0  ♦  4»HS- 

2  *  F  ES2 ( SOL  >  -  2*H+  ♦  2*E 

27 

1 

| 

1 

E  = 

0.1481 

-  0  0507PH 

1 

j 

| 

FE203  4  5*H20  ♦  2*HS-  - 

2  *  +  12*H*  ♦  1 4  *  E - 

28 

1 

1 

E  = 

-07202 

+  0  .  1  1 8  3PH 

1 

1 

| 

HFE02-  -  H20 

F EOH  +  ♦  -  2  *H+  4  E- 

27 

1 

1 

E  = 

-0.6245 

♦  0 . 0592PH 

1 

1 

HFE02- 

F E  <  OH  >24  -  H+  4  E- 

30 

1 

1 

1 

PH  = 

7.58 

| 

1 

F  E  CL  2 ( AQ  >  4  2  *  H2  O  -  2*CL- 

HFE02-  «■  3  * H  + 

3  X 

1 

1 

1 

E- 

-1  .0730 

♦  0  .  1  7  75PH 

1 

1 

HFE02-  -  2  *H20  +  CL-  « 

F  E  CL  ♦  -  3  *  H*  ♦  E- 

32 

t 

1 

1 

PH  = 

14.81 

1 

1 

F  E  S  <  SOL )  ♦  2  *H20  -  HS-  « 

HFE02-  4  2  *H4 

33 

1 

1 

i 

E  = 

-  1  .25  38 

4  0.0276  PH 

1 

1 

HFE02-  -  2  *H20  4  2*HS- 

FES2 ( SOL )  -  H*  4  2  *  E- 

34 

I 

i 

! 

E  = 

0  0270 

-  0 . 0444PH 

1 

i 

HFE02-  4  2  *H20  4  HS- 

*  6  * H+  ♦  8  *  E- 

35 

i 

1 

i 

PH  = 

5.00 

; 

i 

FEOH44  *  H20  = 

F  E  <  OH )  2  ♦  H  + 

36 

1 

1 

i 

E- 

0.7776 

-  0.0572  PH 

i 

i 

FECL2 ( AQ )  4  H20  -  2«CL- 

FEOH44  ♦  Hv  4  E- 

37 

i 

i 

i 

PH- 

2.58 

i 

i 

FECL44  ♦  H20  -  CL- 

FEOH+t  4  H+ 

38 

i 

i 

i 

Ci 

0.8322 

i 

! 

FES(SOL)  4  H20  -  HS- 

FEOH44  4-  E- 

37 

i 

i 

i 

E- 

-1  .5874 

-  0.05 7 2  PH 

i 

i 

FEOH44  -  H20  4-  2  *HS- 

FES2  <  SOL  >  4-  H4-  4-  E- 

40 

i 

i 

i 

E- 

0.1646 

-  0.0676  PH 

i 

i 

F EOH ♦  ♦  ♦  3*H20  4-  HS-  - 

♦  8*H4-  ♦  7  *  E- 

41 

i 

i 

! 

E  — 

1.0753 

-  0  .  1  1 83PH 

i 

i 

F  ECL  2 ( AQ )  4  2  *H20  -  2«CL-  - 

F E ( OH  >24  4  2*H4  4  E- 

42 

1 

1 

i 

PH- 

3.77 

i 

j 

FECL44  4  2 *H20  -  CL-  - 

FE ( OH  >  24  4  2«H4 

43 

{ 

i 

i 

E- 

l  .  127? 

-  0  . 0  5  7  2  PH 

i 

i 

FES(SOL)  4  2  *H20  -  HS-  - 

F  E  <  OH  >24  4  H4  4  E- 

44 

i 

i 

i 

E- 

-1.8831 

i 

i 

i 

F E ( OH ) 2  4  -  2 *H20  4  2«HS-  - 

F  ES2 ( SOL  >  4  E- 

4  S 

i 

! 

i 

E- 

0.1223 

-  0 .0592  PH 

i 

i 

i 

f 

i 

i 

F E ( OK >  24  4  2 *H20  4  HS- 

4  7*H4  4  7  *  E- 

46 

i 

i 

j 

E« 

0 . 6268 

FECL2 ( AQ )  -  CL- 

FECL44  4  E- 

I 


Mb  - 


<1  7 

PH  = 

CD 

O 

1 

1 

| 

1 

ftCL'ZlAUI  -  Z"l-L-  ♦  Mb-  » 

FES<  SOL )  ♦  H* 

4  8 

E  = 

-  0  4  0  3  9 

-  0  0  5  9  2  PH 

1 

1 

| 

F  E  C  L  2  <  AQ  >  -  2  *  CL  -  ♦  2«HS- 

F  E  S  2 ( SOL )  ♦  2  *  H*  ♦  2«E- 

4  9 

E  = 

0  2  4  1  4 

-  0  0666  PH 

1 

1 

| 

F  E  C  L  2 ( AQ )  ♦  4  *  H20  -  2*CL-  ♦  HS 

+  9  *  H  ♦  +  .  8  *  E- 

50 

E  = 

0.6  794 

+  0  0  5  9  2  PH 

1 

1 

| 

FES(SOL)  ♦  CL-  -  HS- 

F  ECL  +  *  -  H+  «•  E- 

5  1 

E  = 

-  1  .4  3  4  6 

-  0  1  1 8  3  PH 

1 

1 

| 

F  ECL ♦  +  -  CL-  ♦  2  *  HS- 

FES2(SOL>  *  2  *  H+-  E- 

52 

E  = 

0.1864 

-  0  .  0  7  6  1  PH 

1 

j 

| 

F  E CL ♦ ♦  ♦  4  *  H2 O  -  CL-  ♦  HS- 

♦  9»H+  ♦  7  *  E - 

5  3 

E  = 

-0.3776 

-  0  . 0  2  9  6  PH 

1 

1 

) 

FES(SOL)  ♦  HS- 

FES2 ( SOL)  +  H*  ♦  2»E- 

54 

E  = 

0.2480 

-  0  .  0  5  9  2  PH 

I 

1 

| 

FES(SOL)  ♦  4  *  H2  O 

♦  8  *  H  ♦  ♦  8  *  E  - 

55 

E  = 

0.4566 

-  0 . 0690PH 

I 

i 

1 

FES2<SOL>  ♦  4  *  H  20  -  HS- 

♦  7  *  H  ♦  ♦  6  *  E  - 

«  AREA  OF  PREDOMINANCE  FOR:  FE  * 

************************************************ 

NO  LOWER  LIMIT 


!  POTENTIAL  1  !  PH  l  !  POTENTIAL  2  !  PH  2  !  LINE  NO.  ! 


!  -0.527 

«  -0.527 

!  -0.991 

MINIMUM  PH 
MAXIMUM  PH 


!  -2.00  < 
»  -0.89  ! 

!  14  .  8  I  ! 

-2.00 
16.00 


-0.527  !  -0.89! 
-0991  !  14.81  • 
-1.097  !  16.00  ! 


6  ! 

8  ! 

3  ! 


*  AREA  OF  PREDOMINANCE  FOR:  FE304 


!  POTENTIAL  1  !  PH  1  !  POTENTIAL  2  !  PH  2  !  LINE  NO.  ! 


-0.801 

!  14.30 

-0.750 

!  16.00  ! 

18 

1 

-0.801 

!  14.30  ! 

-0.876 

'.  14.81  ! 

1  7 

! 

-0.876 

!  14.81  ! 

-0.841 

!  16.00  ! 

12 

i 

-0.750 

<  16.00  ! 

-0.841 

!  16.00  ! 

i 

*  AREA  OF  PREDOMINANCE  FOR.  FE203  * 


NO 


AREA 


OF 


PREDOMINANCE 


AREA  OF  PREDOMINANCE  FOR:  HFE02- 


! 

POTENTIAL  1 

!  PH  1 

POTENTIAL  2  ! 

PH  2  ! 

LINE  NO. 

1 

-0.876 

!  14.81 

-0.841  ! 

1  6 

00  ! 

1  7 

1 

-0.991 

!  14.81 

-1.097  ! 

1  6 

00  ! 

3 

| 

-0.876 

!  14.81 

-0991  ! 

1  4 

8  1  ! 

32 

I 

-0  .  8  4  1 

■  16.00 

-1.097  ! 

1  6 

00  ! 

1 

*  AREA  OF  PREDOMINANCE  FOR:  FEOH++  * 

************************************************ 


NO  AREA  OF  PREDOMINANCE 


«***««******«*«»***«**««*****«*«***«*****»««**«* 

*  AREA  OF  PREDOMINANCE  FOR:  FE<OH>2*  * 

************************************************ 


NO  AREA  OF  PREDOMINANCE 


*  AREA  OF  PREDOMINANCE  FOR:  FECL2CAQ)  * 

************************************************ 


POTENTIAL  1 

!  ?H  1 

POTENTIAL  2  ! 

PH  2 

LINE  NO. 

1 

-0.286 

!  -2.00 

-0.351  ! 

-0.89 

48 

1 

-0.527 

!  -2.00 

-0.527  ! 

-0.89 

6 

1 

-0.286 

!  -2.00 

-0.527  ! 

-2.00 

1 

-0.351 

!  -0.89 

-0.527  ! 

-0.89 

47 

1 

*  AREA  OF  PREDOMINANCE  FOR:  FECL+*  * 


NO  AREA  Of  PREDOMINANCE 


AREA  OF  PREDOMINANCE  FOR  F  E  S ( SQL  >  * 


1 

POTENTIAL  1  ! 

PH 

1 

POTENTIAL 

2  ! 

PH  2 

1 

LINE  NO 

1 

1 

-0.351  ! 

-0 

89 

-0.801 

1 

1  4 

30 

I 

53 

1 

! 

-0.801  ! 

1  4 

30 

-0.876 

1 

1  4 

8  1 

1 

1  7 

1 

! 

-0.527  ! 

-0 

8  9 

-0.991 

I 

14.81 

1 

8 

| 

! 

-0.351  ! 

-0 

8  9 

-0.527 

j 

-0 

89 

1 

4  7 

1 

j 

-0.878  ! 

1  4 

8  1 

-0 . 991 

| 

1  4 

8  1 

1 

32 

1 

«  AREA  OF  PREDOMINANCE  FOR:  FES2<SOL>  * 

************************************************ 


1 

POTENTIAL  1  ! 

PH  1 

POTENTIAL  2  ! 

PH  2 

LINE  NO. 

0.595  ! 

-2.00 

-0.848  ! 

18.00 

55 

-0.288  ! 

-2.00 

-0.351  ! 

-0.89 

48 

-0.351  ! 

-0.89 

-0.801  ! 

14.30 

53 

-0.801  ! 

14.30 

-0.750  ! 

18.00 

1  8 

0.595  ! 

-2.00 

-0.285  ! 

-2.00 

-0.848  ! 

18.00 

-0.750  ! 

18.00 

a*********************************************** 

*  AREA  OF  PREDOMINANCE  FOR:  * 

****•*•******•*••*************•••*•*•••••*••**** 


NO  UPPER  LIMIT 


I  POTENTIAL  1  !  PH  1  I  POTENTIAL  2  !  PH  2  !  LINE  NO.  ! 


!  0 . 5?5 

MINIMUM  PH  - 
MAXIMUM  PH  - 


!  -2.00!  -0648 

-2.00 
18.00 


18.00!  55  ! 


AREA  OF  PREDOMINANCE  FOR:  HFE02- 


•  DISSOLVED  SPECIES  DIAGRAM 


NO  LOWER  LIMIT 


1  POTENTIAL  t  !  PH  1  ‘  POTENTIAL  2  '  PH  2  '  LINE  NO  ' 

!  -0  396  !  9  58  !  -0.681  1  16  00  '  34  ! 

MINIMUM  PH  »  9.58  LINE  NO.  30 

MAX l MUM  PH  «  16.00 


*  AREA  OF  PREDOMINANCE  FOR.  FEOH++  *  DISSOLVED  SPECIES  DIAGRAM 


NO  AREA  OF  PREDOMINANCE 


*  AREA  OF  PREDOMINANCE  FOR:  FE(OH>2«-  *  DISSOLVED  SPECIES  DIAGRAM 


NO  AREA  OF  PREDOMINANCE 


AREA  OF  PREDOMINANCE  FOR:  FECL2 ( AQ  >  *  DISSOLVED  SPECIES  DIACRAM 


NO  LOWER  LIMIT 


•  POTENTIAL  1  !  PH  l  !  POTENTIAL  2  !  PH  2  !  LINE  NO. 

!  0.373  !  -2.00  !  -0.396  !  9.58  !  49 

MINIMUM  PH  =  -2  .  00 

MAXIMUM  PH  =>  9.58  LINE  NO.  30 


*  AREA  OF  PREDOMINANCE  FOR:  FECL++  *  DISSOLVED  SPECIES  DIACRAM 


NO  AREA  OF  PREDOMINANCE 


«  AREA  OF  PREDOMINANCE  FOR:  «  DISSOLVED  SPECIES  DIACRAM 


NO  UPPER  LIMIT 


I  POTENTIAL  II  PH  1  I  POTENTIAL  21  PH  2  1  LINE  NO.  ! 

I  0.373  1  -2.00  !  -0.396  I  9.58  I  4*  I 

I  -0.396  I  9.38  I  -0.681  I  16.00  I  34  ! 

MINIMUM  PH  -  -2.00 

MAXIMUM  PH  «  16.00 


